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Preface 


The  aim  of  this  experiment  was  to  determine  the  mobility  of  an  ion  in  a  buffer  gas, 
specifically  N2^  in  helium,  via  spectroscopic  means.  That  this  experiment  failed  to  produce 
a  value  for  the  mobility,  due  to  the  insensitivity  of  a  the  fitting  function  to  the  parameter  of 
interest,  does  not  imply  a  failure  of  the  experiment.  The  lessons  learned  from  this 
experiment  underscore  the  necessity  of  understanding  the  statistical  significance  of 
experimental  results,  and  the  necessity  of  using  statistical  methods  to  design  an 
experiment.  These  lessons  are  particularly  timely  as  I  prepare  to  assume  my  duties  in  a 
gas  analysis  laboratory,  where  I  will  be  involved  in  improving  the  performance  of  a  gas 
analysis  system. 

I  would  like  to  thank  Dr.  Charles  DeJoseph,  my  advisor  at  Wright  Laboratory,  for 
allowing  me  to  work  with  his  beloved  Bomem  Fourier  Transform  Spectrometer,  and 
allowing  me  to  team  by  example  during  the  summer  months.  Thanks  also  go  to  Dr. 
William  Bailey,  my  AFIT  advisor,  for  influencing  my  selection  of  this  project,  and  helping 
me  achieve  closure  with  this  thesis  by  assisting  me  with  his  insight  of  the  kinetics  of  the 
discharge.  I  am  also  grateful  to  Dr.  Alan  Garscadden,  of  Wright  Laboratory,  for  serving 
on  the  committee,  and  for  offering  possible  explanations  for  the  disagreement  between  my 
measurements,  and  those  of  other  researchers.  (Alan,  after  five  years,  you  finally  got  a 
presentation  out  of  me!).  I  am  indebted  to  Dr.  Rajesh  Nagpal,  a  participant  in  the 
AFOSR  Senior  Investigator  Program,  who  provided  the  calculations  of  the  electron 
energy  distribution  fimction,  which  were  essential  for  explaining  the  discrepancy  between 
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my  observations  and  previously  published  results.  I  am  also  indebted  to  Dr.  Kirk  Mathews 
for  providing  guidance  in  solving  the  nonlinear  least-squares  fitting  problem.  I  would  also 
like  to  recognize  the  significant  contributions  of  Mike  Ray,  for  crafting  four  different 
discharge  tubes  for  this  study,  and  Robert  Knight,  for  providing  me  with  the  use  of  his 
pumpstation. 

On  a  personal  note,  I  would  like  to  take  a  moment  to  thank  those  who  contributed 
to  my  scientific  career.  Success  in  the  pursuit  of  excellence  is  the  result  of  responding  to 
the  challenge  of  high  standards.  Some  who  set  high  standards  in  my  life  were  not 
appreciated  at  the  time,  so  now  I  set  the  record  straight:  Sr.  Lucinda  Girard,  my  high 
school  algebra  instructor,  deserves  special  thanks  for  laying  the  foundation  of  my 
knowledge  of  mathematics.  Her  method  was  simple  and  effective:  homework, 
homework,  homework,  and  when  you’re  tired  of  it,  more  homework.  Mr.  Stephen 
Hamersky  and  Bro.  Timothy  Pieprzyca  also  placed  high  demands  in  their  physics  and 
chemistry  classes,  and  strongly  influenced  my  choice  of  career.  However,  none  of  these 
opportunities  would  have  been  available,  were  it  not  for  the  love,  inspiration,  and  financial 
sacrifice  of  my  parents,  Eugene  and  Mary  Kee,  who  were  determined  to  provide  a  quality, 
value-based  education  for  their  six  children.  To  them,  this  work  is  dedicated. 


Patrick  D.  Kee 


Table  of  Contents 


Page 


Preface . i 

List  of  Figures . v 

List  of  Tables . viii 

Abstract . ix 

I.  Introduction . 1 

A.  Purpose  and  Relevance . 1 

B.  Definition  of  Ion  Mobility . 1 

C.  Approaches  to  the  Experimental  Determination  of  Ion  Mobility . 2 

1 .  Drift  Tube  Measurements . 2 

Biondi  and  Chanin . 3 

McFarland,  et  al . 4 

2.  Optical  Measurements  of  Drift  Velocity . 4 

Haese,  Pan,  and  Oka . 5 

Gudeman,  et  al . 6 

Radunsky  and  Saykally . 7 

Martin  and  Guelachvili . 7 

Hong  and  Miller . 7 

D.  Problem,  Scope,  and  Approach . 11 

E.  Overview . 12 

U.  Theory . 13 

A.  Determination  of  Mobility . 13 

1.  Calculation  of  Time  Averaged  Observed  Ensemble  Velocity . 14 

2.  Derivation  of  Ion  Average  Velocity . 1 5 

3 .  Derivation  of  Emission  Probability  Density . 15 

4.  Derivation  of  E-scaled  Doppler  shift . 17 

5.  Calculation  of  Ion  Mobility  from  Collision  Cross-Section . 17 

B.  Theory  of  Fourier  Transform  Spectrometer . 19 

1 .  The  Interferogram  and  its  Relation  to  the  Spectrum . 19 

2.  Sampling  Function . 21 

3 .  Recovery  of  Spectra . ; 22 


ui 


4.  Elimination  of  Phase  Error . 23 

5.  Limits  of  Resolution . 23 

III.  Equipment  and  Software . 25 

A.  Discharge  Tube . 25 

B.  Vacuum  System  (Including  Flow  System) . 27 

1.  Static  Gas  Fill . 27 

2.  Flowing  Gas  System . 28 

C.  DC  Discharge  Power  Supply . 29 

D.  Fourier-Transform  Spectrometer . 29 

E.  Electric  Field  Measurement  System . 32 

IV.  Procedure . 34 

A.  Measurement  of  Tube  Dimensions . 34 

B.  Calibration  of  FTS  System:  Phase  File . 34 

C.  Alignment  of  Discharge  Tube . 35 

D .  Acquisition  of  Spectra . 36 

E.  Determination  of  Ion  Temperatures . 38 

F.  Voltage  Measurements . 39 

V.  Results  and  Analysis . 40 

A.  Data . 40 

1.  Static  Gas  Fill . 40 

Electric  Field . 40 

Doppler  Shift . 42 

Temperature . 44 

2.  Flowing  Gas . 47 

Electric  Field . 47 

Doppler  Shift . 49 

Temperature . 50 

B.  Analysis . 54 

1 .  Corrections  to  Electric  Field  in  Low  Pressure  Limit . 54 

2.  E-scaled  Doppler  Shift  vs  Pressure . 57 

3 .  Disagreement  with  Hong-Miller  Experiment . 61 

4.  Plot  of  surface . 63 

Comparison  to  Hong-Miller  results . 71 

VI.  Conclusion  and  Recommendations . 75 

A.  Summary  of  Results . 75 

B.  Recommendations  for  Further  Research . 75 


IV 


Appendix  A;  Example  Search  for  Minima  of  . 77 

Appendix  B:  Residence  Time  in  Discharge  Tube . 81 

Appendix  C;  Resolution  of  Spectrometer . 83 

Bibliography . 87 

Vita . 90 


List  of  Figures 


Figure  Page 


1 .  Biondi-Chanin  Experiment . 3 

2.  Haese,  Pan,  and  Oka  Experiment . 6 

3 .  Hong  and  Miller  Experiment . 8 

4.  Doppler  Shift  in  N2^  emission  spectra .  9 

5.  Michelson  Interferometer . 20 

6.  Discharge  Tube . 25 

7.  Probe  Design . 26 

8.  Static  Gas  Fill  Vacuum  System . 28 

9.  Schematic  of  Fourier  Transform  Spectrometer . 31 

10.  Electric  Field  Measurement  System . 32 

1 1 .  Peak  Picking  Program . 33 

12.  E  vs  P,  Static  Fill,  I  =  1  mA . 40 

13.  E  vs  P,  Static  Fill,  1  =  5  mA . 41 

14.  E  vs  P,  Static  Fill,  1=10  mA . 42 

15.  Doppler  Shift  vs  Pressure,  Static  Fill,  1=1  mA . 43 

16.  Doppler  Shift  vs  Pressure,  Static  Fill,  1  =  5  mA . 43 

17.  Doppler  Shift  vs  Pressure,  Static  Fill,  1=10  mA . 44 

18.  Ion  Temperature  vs  Pressure,  Static  Fill,  I  =  1  mA . 45 

19.  Ion  Temperature  vs  Pressure,  Static  Fill,  1  =  5  mA . 45 

20.  Ion  Temperature  vs  Pressure,  Static  Fill,  1=10  mA . 46 

2 1 .  Ion  Temperature  vs  Current,  Static  Fill . 46 

22.  E  vs  P,  Flowing  Gas,  1=1  mA . 47 

23 .  E  vs  P,  Flowing  Gas,  1  =  5  mA . 48 

24.  E  vs  P,  Flowing  Gas,  1=10  mA . 48 

25.  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1=1  mA . 49 

26.  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1  =  5  mA . 50 

27.  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1=10  mA . 50 

28.  Ion  Temperature  vs  Pressure,  Flowing  Gas,  I  =  1  mA . 51 

29.  Ion  Temperature  vs  Pressure,  Flowing  Gas,  1  =  5  mA . 52 

30.  Ion  Temperature  vs  Pressure,  Flowing  Gas,  1=10  mA . 52 

3 1 .  Ion  T emperature  vs  Current,  Static  Fill . 53 

32.  Adjusted  Electric  Field  vs  Pressure,  Static  Fill,  1=1  mA . 54 

33 .  Adjusted  Electric  Field  vs  Pressure,  Static  Fill,  1  =  5  mA . 55 

34.  Adjusted  Electric  Field  vs  Pressure,  Static  Fill,  1=10  mA . 55 

35.  Adjusted  Electric  Field  vs  Pressure,  Flowing  Gas,  1=1  mA . 56 

36.  Adjusted  Electric  Field  vs  Pressure,  Flowing  Gas,  1  =  5  mA . 56 

37.  Adjusted  Electric  Field  vs  Pressure,  Flowing  Gas,  1=10  mA . 57 


VI 


38.  E-scaled  Doppler  Shift  vs  Pressure,  Static  Fill,  1=1  mA . 58 

39.  E-scaled  Doppler  Shift  vs  Pressure,  Static  Fill,  1  =  5  mA . 58 

40.  E-scaled  Doppler  Shift  vs  Pressure,  Static  Fill,  1=10  mA . 59 

41 .  E-scaled  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1=1  mA . 59 

42.  E-scaled  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1  =  5  mA . 60 

43.  E-scaled  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1=10  mA . 60 

44 .  N2"^(B)  F  ormation  Processes . 61 

45 .  Electron  Energy  Distribution  Function . 62 

46.  Plots  Static  Fill,  1=1  mA . 66 

47.  ^  Plot,  Static  Fill,  1  =  5  mA . 67 

48.  Plot,  Static  Fill,  1  =  10  mA . 68 

49.  x^  Plot,  Flowing  Gas,  1  =  1  mA . 69 

50.  x^  Plot,  Flowing  Gas,  1  =  5  mA . 70 

51.  x^  Plot,  Flowing  Gas,  1=10  mA . 71 

52.  x^  Plot,  Hong-Miller  Experiment,  1  =  25  mA . 73 

53.  Insensitivity  of  Fitting  Equation  to  Choice  of  Qv  for  Hong-Miller  Data . 74 


vii 


List  of  Tables 


Table  Page 


1.  values.  Static  Fill,  I  =  1  mA . 66 

2.  values.  Static  Fill,  1-5  mA . 67 

3.  x^  values.  Static  Fill,  1=10  mA . 68 

4.  x^  values.  Flowing  Gas,  1=1  mA . 69 

5.  x^  values.  Flowing  Gas,  1  =  5  mA . 70 

6.  x^  values,  Flowing  Gas,  1=10  mA . 71 

7.  x^  values,  Hong-Miller  Experiment,  I  =  25  mA . 73 


viii 


AFIT/GAP/EN/94D-4 


Abstract 

Recently,  Hong  and  Miller  used  Fourier  Transform  Emission  (FTE)  spectroscopy 
to  determine  the  N2^  mobility  in  helium  by  relating  the  observed  Doppler  shift  of  the  N2^ 
(B-X)  transition  to  the  measured  electric  field  in  the  discharge  as  determined  by  probes. 
In  order  to  establish  the  technique  in  house,  this  effort  duplicated  that  experiment  and 
investigated  the  sensitivity  and  utility  of  that  method.  Doppler  shift  measurements  in  this 
experiment  disagreed  significantly  with  those  reported  by  Hong  and  Miller.  Attempts  to 
reduce  the  data  (and  thus  determine  the  mobility)  using  a  model  developed  by  Hong  and 
Miller  were  unsuccessful.  The  model  did  not  qualitatively  or  quantitatively  describe  the 
Doppler  shifts  observed  in  this  work.  At  this  time,  an  explanation  for  the  differences 
between  the  two  experiments  is  not  known.  Possibilities  include  effects  of  striations  on 
the  Doppler  shift  measurements,  competition  between  one-  and  two-step  excitation  of  the 
N2^  (B)  state,  and  failure  to  interpret  the  Hong  and  Miller  model  correctly.  Subsequent 
attempts  to  model  Hong  and  Miller’s  own  data  (taken  from  their  figures)  yielded  mobility 
parameters  which  disagreed  with  their  published  results. 
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MEASUREMENT  OF  N2^  MOBILITY  IN  He  VIA  FOURIER 
TRANSFORM  EMISSION  SPECTROSCOPY 

I.  Introduction 

A.  Purpose  and  Relevance 

The  study  of  the  transport  of  ions  in  gases  yields  information  applicable  to 
scientific  investigations  of  interest  to  the  Air  Force.  The  determination  of  the  ion  mobility, 
a  coefficient  that  relates  an  ion’s  drift  velocity  in  a  buffer  gas  to  an  applied  electric  field, 
can  yield  collision  cross-section  data  essential  for  a  quantitative  understanding  of  gas 
discharges  and  atmospheric  phenomena  [1:8],  Current  Air  Force  research  projects  that 
may  benefit  from  ion-mobility  data  include;  studies  of  atmospheric  phenomena,  and  gas 
discharge  lasers. 

In  this  effort,  the  mobility  of  N2^  ions  will  be  measured  in  a  carrier  gas  using 
Fourier-transform  spectroscopy,  and  will  be  compared  to  previously  published  values 
obtained  from  both  drift  tube  and  spectroscopic  measurements. 

B.  Definition  of  Ion  Mobility 

The  ion  mobility,  K,  is  defined  as: 


1 


where  va  is  the  drift  velocity  of  an  ion  in  either  a  parent  or  buffer  gas,  and  E  is  the  electric 
field.  The  measurements  of  the  two  quantities  are  necessary  for  the  determination  of  the 
ion  mobility.  Electric  field  measurements  are  often  made  using  probes,  while  drift  velocity 
measurements  have  been  made  by  several  methods.  These  measurements  can  be 
conveniently  divided  into  two  areas;  drift  tube  measurements,  and  optical  drift  velocity 
measurements.  Illustrative  examples  of  both  techniques  are  discussed  in  the  following 
sections. 

Ion  mobility  measurements  are  taken  under  a  variety  of  experimental  conditions. 
To  compare  the  ion  mobilities  taken  under  different  temperatures  and  pressures,  a  quantity 
known  as  the  reduced  ion  mobility,  Ko,  is  defined  as: 


where  K,  P,  and  T  are  respectively  the  experimental  values  for  the  ion  mobility,  buffer  gas 


pressure,  and  buffer  gas  temperature. 


C.  Approaches  to  the  Experimental  Determination  of  Ion  Mobility 

1 .  Drift  Tube  Measurements.  A  drift  tube  has  been  the  most  common  device 
used  to  measure  the  drift  velocity  of  ions.  The  principle  of  operation  is  to  measure  the 
time  of  flight  of  ions  through  a  buffer  gas,  across  a  known  distance.  The  essential 
characteristics  of  a  drift  tube  apparatus  are  a  vacuum  system,  a  pulsed  ion  source,  a  pair  of 
electrodes  to  accelerate  the  ion  cloud,  probes  for  measuring  the  electric  field,  and  an 
ammeter  connected  to  the  collector  [1:32],  Variations  in  drift  tube  design  may  include 
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mass  spectrometers  to  separate  ions  of  interest  from  other  ionic  species,  and  flowing- 
afterglow  discharges  for  ion  production  for  the  drift  tube. 

Biondi  and  Chanin.  In  an  experiment  by  M.A.  Biondi  and  L.M.  Chanin 
[2],  (see  Figure  1),  ions  are  produced  by  a  high  voltage  pulse  in  a  drift  tube  in  a  region 
between  an  electrode  and  a  grounded  grid.  An  electric  field  applied  between  the  electrode 
and  the  grid  accelerates  the  ions,  which  penetrate  the  grid  and  enter  the  drift  space  of  the 
tube.  The  ion  current  in  the  drift  tube  induces  a  current  across  a  resistor  in  an  external 
circuit.  The  voltage  across  the  resistor  is  measured  and  displayed  on  an  oscilloscope, 
with  an  external  trigger  supplied  by  the  pulse  generator.  Then  with  va  =  d/t  and  E  =  V/d, 
the  ion  mobility  is  determined  by. 

K  =  —  (3) 

Vt 

where  d  is  the  spacing  between  the  grid  and  the  negative  electrode,  V  is  the  potential 
difference  between  the  grounded  grid  and  the  negative  electrode,  and  t  is  the  transit  time 
of  the  ion  across  the  drift  space.  Similar  techniques  were  used  to  determine  the  N2^  ion 
mobility  in  helium  [3],  yielding  a  result  of  K«  =  19  ±  2  cm^/volt/sec. 


Discharge  Region 

_ ^ 

Drift  Region 

(- 

) 

1 

Grid 

- > 

E  field 

Collector  Electrode 

Figure  1.  Biondi-Chanin  Experiment 


3 


McFarland,  et  al.  McFarland,  et  al,  combined  a  flowing  afterglow  system 
with  a  drift  tube  apparatus  to  measure  the  mobilities  of  a  variety  of  ions,  including  N2^, 
with  helium  as  the  buffer  gas  [4],  The  flowing  afterglow  system  was  used  for 
measurement  of  reaction  rates  by  accelerating  a  steady-state  stream  of  ions  through  a 
reactant  gas  in  a  drift  tube.  When  the  reactant  gas  was  replaced  with  a  buffer  gas,  and  the 
ions  produced  in  a  pulsed  mode,  the  ion  mobility  was  then  determined  from  the  time  of 
flight  of  the  ion  through  the  drift  tube.  The  results  of  their  measurement  yielded  a  value  of 
Ko  =  21.0±  1.3  cm^/volt/sec,  which  agrees  with  the  previously  measured  value  of 

Johnsen,  Brown,  and  Biondi  [3].  Advantages  of  the  flowing  afterglow/drift  tube 
apparatus  over  a  conventional  drift  tube  device  include  the  ability  to  measure  the  mobility 
of  negative  ions  (which  was  done  with  H*,  O',  and  OH'  in  helium).  In  addition, 
measurements  were  made  of  the  reaction  rates  of  ions  with  neutrals  as  a  function  of  the 
ion  kinetic  energy. 

2.  Optical  Measurements  of  Drift  Velocity.  Advances  in  spectroscopy  in  the  last 
fifteen  years  have  permitted  the  observation  of  high-resolution  vibration-rotation  spectra 
of  gas  phase  molecular  ions  in  a  dc  discharge  [5].  If  a  spectrometer  has  sufficient 
resolution  to  measure  the  spectral  shift  due  to  the  Doppler  effect  of  an  ion  in  a  dc 
discharge,  then  the  drift  velocity  of  the  ion  can  be  determined  from  the  Doppler  shift.  The 
advantage  of  this  technique  over  a  drift  tube  apparatus  is  that  the  spectra  of  an  ionic 
species  can  be  distinguished  from  other  ions  in  the  discharge,  whereas  with  a  drift  tube, 
care  must  be  taken  to  differentiate  among  the  time-of-arrival  of  many  ionic  species. 
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The  first  successful  spectroscopic  measurement  of  an  ion  mobility  was  performed 
by  Haese,  Pan,  and  Oka  using  infrared  laser  spectroscopy  [5],  Subsequent  work  in 
detecting  faint  molecular  ion  absorption  lines  was  pioneered  by  Gudeman,  et  al,  with  a 
technique  known  as  ‘velocity  modulation’  [6],  Radunsky  and  Saykally  then  applied 
velocity  modulation  to  the  measurement  of  axial  electric  fields  in  a  He/N2  discharge  by 
using  the  established  value  of  the  N2^  mobility  in  helium  and  the  observed  Doppler  shift  in 
the  N2^  spectra  [7],  Velocity  modulation  was  extended  to  the  case  of  emission  by  Martin 
and  Guelachvili  [8],  who  combined  the  technique  with  Fourier-transform  spectroscopy 
(FTS)  for  the  selective  detection  of  molecular  ions,  and  discussed  the  feasibility  of 
applying  this  method  to  the  measurement  of  quantum-state  specific  ion  mobilities.  Most 
recently,  Hong  and  Miller  applied  Fourier-transform  spectroscopy  with  a  reversible- 
polarity  discharge  with  results  that  agree  well  with  previous  results  [9], 

Haese.  Pan,  and  Oka.  Haese,  Pan,  and  Oka  performed  the  first  successful 
spectroscopic  ion  mobility  measurement,  using  infrared  laser  spectroscopy  to  observe  the 
drift- velocity  induced  Doppler  shift  in  ArH^  known  vibration-rotation  absorption  lines  [5]. 
Using  a  beamsplitter,  Haese,  Pan,  and  Oka  illuminated  the  positive  column  of  the 
discharge  from  opposite  ends  of  the  discharge  tube,  simultaneously  measuring  the  blue  and 
red  shifted  absorption  spectra  (see  Figure  2).  The  drift  velocity  of  the  ions  was  deduced 
from  the  Doppler  shift  in  the  vibration-rotation  absorption  lines,  while  the  axial  electric 
field  was  measured  with  platinum  wire  probes.  The  calculated  values  for  the  ion  mobility 
of  ArH"  agreed  with  a  previous  mass-spectroscopic  drift-tube  study  [10]. 
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Figure  2.  Haese,  Pan,  and  Oka  Experiment 

Gudeman,  et  al.  One  complication  faced  in  obtaining  ion  spectra  when  the 
ion  number  density  is  a  small  fraction  of  the  total  gas  density,  is  that  the  signal  from  ion 
absorption  may  be  much  less  than  that  from  neutral  absorption  [6],  One  of  the  first 
successful  measurements  of  the  infrared  spectrum  of  a  molecular  ion  was  by  Gudeman, 
Begemann,  Pfaff,  and  Saykally,  using  the  “velocity  modulation”  method  [6],  The  HCO^ 
ion  absorption  spectrum  was  obtained  by  modulating  the  polarity  of  a  gas  discharge  at 
several  kilohertz.  The  modulation  of  the  discharge  polarity  resulted  in  both  red-  and  blue- 
shifted  absorption  peaks.  This  Doppler  shift  was  used  to  separate  ion  transitions  from  the 
more  intense  neutral  transitions.  The  ion  mobility  was  not  measured  in  this  experiment, 
but  the  concept  of  modulating  the  polarity  of  the  discharge  to  obtain  ionic  spectra  would 
be  adopted  by  Hong  and  Miller  for  their  ion  mobility  measurements. 
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Radunskv  and  Saykally.  Radunsky  and  Saykally  applied  the  velocity 
modulation  method  in  combination  with  laser  absorption  spectroscopy  to  determine  the 
axial  electric  field  in  a  glow  discharge  of  a  N2/He  mixture  [7],  In  their  experiment,  the 
shift  in  the  absorption  spectra  of  was  measured,  and  related  to  the  axial  electric  field 
using  the  previously  published  value  of  the  N2^  mobility  in  helium. 

Martin  and  Guelachvili.  This  velocity  modulation  method  was  extended 
to  the  case  of  emission  by  Martin  and  Guelachvili  [8],  They  combined  velocity  modulation 
of  the  ions  with  Fourier  transform  spectroscopy  (FTS)  to  selectively  detect  molecular  ions 
in  a  discharge,  and  to  measure  the  ratio  of  the  dipole  moment  to  the  vibrational  transition 
moment  of  ArH"^.  FTS  has  the  advantage  over  laser  absorption  spectroscopy  of  having 
wide  and  continuous  spectral  coverage,  as  well  as  having  more  consistent  wave-number 
accuracy.  No  ion  mobility  measurements  were  made  in  this  experiment,  however,  it  was 
suggested  in  this  paper  that  this  technique  could  be  used  to  determine  the  excited-state 
mobilities  of  ions  [8]. 

Hong  and  Miller.  Hong  and  Miller  applied  the  FTS  technique  to 
determine  the  ion  mobility  of  N2^  in  He,  and  fit  their  data  to  a  theoretical  model  to 
determine  the  velocity  changing  and  quenching  cross  sections  for  He-N2^  collisions  [9]. 
N2^  ions  were  produced  in  a  dc  discharge,  the  Doppler  shift  in  both  the  B^Su^(v=0)- 
X^Eg^(v=l),  and  the  A^nu(v=l)-X^lg^(v=0)  transitions  were  observed.  To  select  red-  or 
blue-shifting  of  the  emission  spectra,  the  polarity  of  the  electrodes  could  be  swapped  by  a 
computer  controlled  switch  connected  to  the  power  supply  (see  Figure  3).  The  discharge 
was  viewed  with  a  high-resolution  Fourier  transform  spectrometer,  and  successive  spectra 
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were  taken  by  first  viewing  from  one  polarity,  then  with  the  polarity  reversed.  Because 
the  spectra  were  taken  under  nearly  identical  conditions  (except  direction  of  ion  drift 
velocity),  the  only  differences  in  the  resulting  spectra  were  the  opposite  signs  of  Doppler 
shift.  Taking  the  difference  of  the  red-  and  blue-shifted  spectra  isolated  the  emission 
peaks  due  to  ionic  transitions,  and  the  Doppler  shift  was  obtained  by  taking  half  the 
separation  between  the  blue  and  red  shifted  peaks.  Figure  4  displays  the  Doppler  shift 
observed  by  the  author  in  the  spectra  as  a  result  of  switching  the  electrode  polarities 
between  measurements. 


Figure  3.  Hong  and  Miller  Experiment. 
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Figure  4.  Doppler  Shift  hi  N2+  emission  spectra  (Author's  data). 


The  electric  field  was  estimated  by  using  one  of  two  methods.  At  pressures  below 
1  torr,  (no  collision  limit),  the  electric  field  can  be  determined  from  Doppler  shift  of  the 
short-lived  B^X  transition  of  N2^,  (for  which  x=66ns)  [11]; 


where  M  is  the  ion  mass,  T  is  the  collision  frequency,  is  the  mean  free  path  for 
quenching  collisions,  V  is  the  ion-molecule  relative  velocity,  x  is  the  radiative  lifetime  of 
the  upper  quantum  state  of  the  ion,  Av  and  Vo  are  respectively  the  shift  in  the  emission 
frequency,  and  the  emission  frequency,  and  c  is  the  velocity  of  light.  In  the  no  collision 
limit,  (N->0),  equation  (4)  becomes; 
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MFC 


(5) 


Eo  = 


MFc  Av 

1  _  ^  ev„ 

1  +  Ft 


If  the  product  Ft  is  small,  then  performing  the  binomial  expansion  on  (1+Fx)'* 


leads  to: 


E 


D 


Me  Av 

e% 


(6) 


At  pressures  below  than  1  torr,  this  is  a  good  approximation  for  the  short-lived 
B->X  transition.  At  pressures  above  1  torr,  this  is  not  a  good  approximation,  nor  is  it  a 
good  approximation  for  the  longer  lived  A-^X  transition. 

Hong  and  Miller  also  used  probes  to  measure  the  electric  field.  The  electric  field 
measured  by  probes  at  low  pressures  was  observed  to  differ  by  a  constant  of  1.5  V/cm 
fi-om  the  result  obtained  from  equation  (6)  using  the  measured  Doppler  shift.  This  offset 
was  attributed  to  sheath  effects  and  probe  positioning.  The  estimated  value  of  the  electric 
field  was  then  taken  to  be  the  probe  measurement  minus  the  offset  value.  This  estimated 
value  was  used  at  pressures  above  1  torr,  and  for  the  A-X  transition,  where  equation  (6) 
no  longer  applies  (due  to  the  longer  lifetime  of  the  A  state).  Hong  and  Miller’s  results  for 
the  electric  field  agreed  well  with  Radunsky  and  Saykally,  who  used  velocity  modulation 
IR  laser  absorption  spectroscopy  in  a  similar  N2-He  glow  discharge  [7].  The  resulting  N2^ 
reduced  ion  mobility  for  the  B-^X  transition  was  16.6  cmW/s,  and  20.7  cm^/volt/sec  for 
the  A^X  transition,  as  compared  to  the  established  value  of  19  cm^V/s  [3]. 
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D.  Problem.  Scope,  and  Approach. 


In  this  effort,  the  mobility  of  the  ion  is  measured  in  helium  using  Fourier- 
transform  emission  spectroscopy.  This  value  for  the  mobility  is  compared  to  results 
obtained  from  both  drift  tube  and  spectroscopic  measurements.  The  N2^  -  He  collision 
cross-section  is  determined  as  a  fitting  parameter  from  a  relationship  between  the  Doppler 
shift  of  the  ionic  spectra  to  the  pressure.  Then  the  mobility  is  computed  from  the  collision 
cross-section.  The  development  of  these  relationships  is  presented  in  the  Theory. 

The  transition  observed  in  computing  the  Doppler  shift  was  the  B^Su^(v=0)- 
X^Sg^(v=0).  This  differs  from  Hong  and  Miller  who  used  the  weaker  v  =  0->l  transition 
[9].  Measurements  were  not  taken  of  the  A-X  transition.  Pressures  in  this  experiment 
varied  from  0.25  torr  to  6  torr,  while  discharge  current  was  varied  from  1  to  10  mA. 
Higher  currents  were  not  used  due  to  observed  sputtering  in  the  tube.  In  comparison,  the 
Hong-Miller  experiment  maintained  the  discharge  at  25  mA. 

The  experimental  setup  is  similar  to  that  of  Hong  and  Miller,  is  described  in 
Chapter  III.  The  discharge  tube  has  the  same  internal  diameter,  however,  a  significant 
variation  is  that  provisions  are  made  in  this  effort  for  both  flowing  gas  and  a  static  gas  fill. 
Axial  electric  fields  are  measured  using  four  bare  tungsten  voltage  probes  inserted  in  the 
positive  column.  The  probe  voltage  is  plotted  as  a  function  of  distance  of  the  probe  from 
an  electrode,  and  the  absolute  value  of  the  slope  is  taken  as  the  electric  field. 
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E.  Overview. 

The  following  sections  present  in  detail  the  effort  to  determine  the  N2^  mobility  in 
helium.  Essential  background  theory  is  developed  and  presented  in  Chapter  II.  This 
includes  the  expressions  that  relate  the  Doppler  shift  in  the  N2^  spectra  to  the  collision 
cross-section  for  N2^  and  He,  and  the  computation  of  the  ion  mobility  from  the  collision 
cross-section.  This  also  includes  a  brief  survey  of  the  principles  of  Fourier-transform 
spectroscopy. 

The  equipment  used  for  the  experiment,  and  the  algorithms  used  in  the  analysis  are 
discussed  in  Chapter  III,  followed  by  a  description  in  Chapter  IV  of  the  experimental  and 
analytical  procedures.  Results  are  analyzed,  presented,  and  discussed  in  Chapter  V. 
Conclusions  drawn  from  this  effort  and  recommendations  for  further  efforts  are  presented 
in  Chapter  VI. 
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II.  Theory 


A.  Determination  of  Mobility 

The  shift  in  the  spectral  lines,  from  an  ensemble  of  ions  drifting  in  a  buffer  gas 
under  the  influence  of  an  electric  field,  can  be  related  to  the  drift  velocity  by  the  Doppler 
relation; 


Av 


(7) 


where  <Ur  >  is  the  time  average  of  the  ensemble  average  velocity  of  the  absorbing/emitting 
ions  along  the  direction  of  the  electric  field.  For  the  case  of  absorption,  the  ground  state 
ion  is  observed,  and  the  Doppler  shift  in  the  absorption  spectra  is  a  measure  of  average 
ground  state  ion  velocity.  This  means  that  <Ur  >  is  equal  to  Vd.  Thus  in  the  case  of 
absorption  spectroscopy,  the  mobility  can  be  directly  calculated  from: 


E  E  v„E 


(8) 


For  emission,  it  is  an  upper  state  that  is  being  observed.  The  probability  that  an 
ion  in  an  excited  state  will  emit  a  photon  is  a  function  of  time,  and  is  related  to  the  lifetime 
of  the  upper  state.  This  finite  lifetime  means  that  the  Doppler  shift  measures  the  time 
average  of  the  mean  velocity  of  the  ions  accumulated  prior  to  emission,  not  the  average 
drift  velocity: 


{«r>  ^  Vd 


(9) 
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Hence,  the  mobility  must  be  calculated  by  another  method  than  equation  (8)  for 
the  case  of  emission  spectroscopy.  To  do  this,  a  simple  model  is  presented  in  the 
following  sections  that  relates  the  observed  Doppler  shift  to  the  collision  cross-section, 
and  thus,  to  the  ion  mobility. 


1.  Calculation  of  Time  Averaged  Observed  Ensemble  Velocity.  The  time  average 
of  the  mean  observed  ion  velocity  for  an  ensemble  of  ions  is  defined  as: 


/  \  io  P(t)Ur(t)dt 

^  '  Io’  P(t)dt 


(10) 


where  P(t)  is  the  probability  of  emission  per  unit  time,  or  emission  probability  density. 
Here,  tq  is  the  mean  time  between  quenching  collisions  for  the  ion  and  is  defined  as: 


t 


q 


(11) 


where  V  is  the  rms  average  velocity  of  the  reduced  collision  mass,  p,  and  A<q  is  the  mean 
tree  path  for  quenching  [9].  The  mean  time  between  quenching  is  chosen  as  the  upper 
limit  for  the  time  average,  because  for  t  >  tq,  ions  are  not  expected  to  be  in  the  excited 
state,  and  no  longer  contribute  to  the  observed  Doppler  shift.  V  is  determined  by  equation 
(12),  where  kg  is  Boltzmann’s  constant,  and  T  is  the  ion/gas  temperature  (assuming 
equilibrium)  [12:121]: 


V 


3kBT 


(12) 
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The  mean  free  path  for  quenching  can  be  expressed  as  [9]; 


^  =  X  Q,(i)Ni  +  ^  (13) 

Aq  1  V 

where  Qq(i)  is  the  quenching  cross  section  for  the  ion  with  the  i^''  buffer  gas  and  yi  is  the 
ambipolar  diffrision  rate  of  arrival  of  ions  at  the  wall. 

2.  Derivation  of  Ion  Average  Velocity.  For  an  ion  moving  under  the  influence  of 
an  electric  fleld  in  a  buffer  gas,  the  differential  equation  of  motion  for  the  ion  can  be 
expressed  as  [13:565]: 

du 

dt  ^  ^  ^ 

where  M  is  the  ion  mass,  aur  is  the  Stoke’ s  resistance  term  and  a  is  a  function  of  the 
viscosity  of  the  buffer  gas. 

Equation  (14)  is  a  linear  first  order  differential  equation,  and  the  solution  is: 

eE 

»r  "  ^  0  “  exp(-rt))  (15) 

where  F  =  a/M  is  the  ion  collision  frequency. 

3,  Derivation  of  Emission  Probability  Density.  The  calculation  of  the  ensemble 
average  velocity  from  equation  (10)  requires  knowing  the  emission  probability  density, 
P(t).  For  emission  from  i->j,  the  probability  density  can  be  expressed  as: 

P(t)  =  pA,,N,  (16) 


IS 


where  Aij  is  the  i— >j  transition  rate,  Ni  is  the  number  of  ions  in  the  excited  state,  and  p  is  a 
proportionality  constant. 

The  number  of  ions  in  the  upper  state  is  a  fiinction  of  time  and  is  given  by 


dN; 

dt 


-Nil 

j 


(17) 


Solving  for  N;  yields: 


Ni  =  N,„  exp( 


(18) 


where 


X 


(19) 


Then  the  probability  density  may  be  expressed  as: 


P(t)  =  exp(-  /^)  (20) 

The  probability  of  emission  over  all  time  must  equal  1,  thus: 

1  =  P(t)dt  =  PA;  exp(-  ^)dt  (21) 


Integrating  equation  (21)  gives: 


X 


(22) 
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And  the  probability  density  is  given  as: 


p(t)  =  -  exp(-  yi) 

T  \  /  X/ 


4,  Derivation  of  E-scaled  Doppler  shift.  Inserting  the  equations  (15)  and  (23) 
into  equation  (10),  and  performing  the  integration  yields  the  time  average  of  the  observed 
mean  ion  velocity; 


1  +  rx 


-  exp  - 


1  -  exp  - 


Dividing  both  sides  of  equation  (24)  by  c  and  E  yields  the  E-scaled  Doppler  shift. 


1  +  Ft 


-  exp  - 


1  -  exp - - 

Vt 


5.  Calculation  of  Ion  Mobility  from  Collision  Cross-Section.  Taking  the  steady 
state  limit  of  equation  (15)  yields  [14]: 


u,  =  V.  = 


The  Wannier  ffee-flight  theory  yields  another  expression  for  the  drift  velocity  for 
an  ion  drifting  in  a  buffer  gas  of  mass  mg  [12:121]; 
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( 


1 


eE 


(27) 


Va  = 


-  +  - 


NQ.V 


where  Qv  is  the  collision  cross  section  for  the  ion/gas  collision,  nig  is  the  mass  of  the  buffer 
gas,  and  ^  is  a  constant  which  is  generally  accepted  to  be  equal  to  1  for  the  case  of  a  heavy 
ion  and  a  light  buffer  gas. 

Combining  equations  (26)  and  (27),  and  solving  for  Y  yields: 


r 


m 


g 


nig  +  M 


NQ,V 


(28) 


This  disagrees  with  Hong  and  Miller,  who  state  for  a  collision  frequency  [9]: 


r 


ni„ 


m^  +  M 


NQ,V 


(29) 


With  a  value  of  Qv  determined  from  fitting  equation  (25)  to  experimental  Doppler 
shift  data  over  a  range  of  pressure,  the  mobility  can  be  determined  from  equations  (26) 
and  (28) 


=  ^  qE  ^  q  M 

E  MEE  MNQ^V  m^ 

(30) 


IlNQ^V 


At  this  point  it  is  worth  mentioning  the  relationship  and  role  of  Qv  and  Qq.  The 
velocity  changing  cross-section,  Qv,  determines  the  steady  state  drift  velocity  of  an  ion. 
The  quenching  cross-section,  Qq,  serves  to  nonradiatively  depopulate  the  upper  state  of 
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the  excited  ion,  and  remove  the  ion  from  the  sample  viewed  by  the  experiment.  Qq  and  Qv 
are  related  in  the  sense  that  a  quenching  collision  is  also  a  velocity  changing  collision.  If 
there  are  nonquenching  velocity  changing  collisions,  then  we  expect  Qv  >  Qq. 

B.  Theory  of  Fourier  Transform  Spectrometer 

The  theory  of  the  Fourier  transform  spectrometer  has  been  well  explained 
elsewhere  [15],  The  following  sections  contain  an  elementary  discussion  of  the 
mathematics  involved  in  the  use  of  a  Fourier  transform  spectrometer  in  this  experiment. 

1  ■  The  Interferogram  and  its  Relation  to  the  Spectrum.  The  essential  components 
of  the  Fourier  transform  spectrometer  are  a  Michel  son  interferometer,  a  detector,  and  a 
computer  for  processing  the  signal  from  the  detector.  Figure  5  is  a  diagram  of  a 
Michelson  interferometer  with  a  detector  at  the  image  plane.  For  this  discussion,  we  will 
assume  an  on-axis,  quasi-monochromatic  point  source  and  a  beamsplitter  of  amplitude 
transmittance  t,  and  a  reflectance  r. 

If  the  incident  wave  is  of  the  form  A  exp[  i(c)t-27rxCT)],  where  ©  and  a  are 
respectively  the  angular  frequency  and  wavenumber  of  the  radiation,  then  the  net 
amplitude  at  the  image  plane  will  be: 

A;,^  =  A(rt)[e‘(“‘-'’^''’>  +  (31) 

where  xi  and  X2  are  respectively  the  roundtrip  distances  from  the  beamsplitter  to  the  fixed 
and  the  movable  mirror. 
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Mirror 


Figure  5.  Michelson  Interferometer 

Then  the  energy  reaching  the  detector  is: 

Edet  =  |Adetf  =  2A>p[l  +  cos(2jt(xl  -  x2)a)]  (32) 

This  consists  of  a  dc  term  and  an  ac  term.  The  ac  term  is  called  the  interferogram. 
Letting  =  B(a)  da,  where  B(a)  is  the  spectral  density,  expressing  I  rtl  ^  as  s,  the 
beamsplitter  efficiency,  and  X1-X2  =  x  as  the  path  difference,  then  the  ac  part  of  equation 
(32)  can  be  expressed  as: 

dl(x)  =  2eB(a)  cos(27rax)  da  (33) 

So  far  we  have  considered  our  source  as  being  quasi-monochromatic.  To 
determine  the  interferogram  resulting  from  a  broad  spectral  range,  we  integrate  dl(x)  with 
respect  to  a: 
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(34) 


I(x)  =  I  dl(x)  =  2ej  B(a)  cos(27iax)da 

o  0 

The  interferogram  is  recorded  as  a  function  of  the  path  difference  x,  by  scanning 
the  movable  mirror  over  its  maximum  path  difference.  The  spectral  density,  B(a)  is  then 
recovered  by  performing  the  inverse  cosine  Fourjer  transform  on  I(x): 


B(a)  =  J  I(x)  cos(27rox)dx 


(35) 


2,  Sampling  Function.  The  signal  received  at  the  detector  is  of  an  analog  nature. 
In  order  to  be  analyzed  by  a  digital  computer,  an  analog-to-digital  conversion  must  take 
place.  This  can  be  represented  by  a  sampling function,  'F  [  1 5] : 


4'(x)  =  I  5(x  -  n) 


where  T'  is  a  series  of  delta  functions  at  regularly  spaced  intervals,  n. 


Some  identities  of  the  sampling  function  are: 


H'(x  +  m)  =  'F(x)  (m  =  any  integer) 


if  n^ 

'P(ax)  =  -  X - 

a  V  a; 


And  the  Fourier  transform  of  is  expressed  as: 


1  ^  (T^ 

FTf'F(ax)]  =  -  'F  - 
a  Va^ 


(36) 


(37) 

(38) 


(39) 
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Then  an  interferogram  which  is  digitally  sampled  may  be  expressed  as; 

l'(x)  =  ■  I(x)  (40) 

VAx; 

And  the  spectrum  obtained  from  performing  a  Fourier  transform  on  the 
interferogram  is: 

'  f  ^ 

B  (a)  =  (Ax)M^  —  *  B(a)  (41) 

VAoj 

where  the  symbol  refers  to  the  convolution  of  the  two  functions,  and  B(ct)  is  the  ‘true’ 
spectrum. 

3.  Recovery  of  Spectra.  In  practical  situations,  interferogram  may  not  be  an  even 
function  of  the  path  difference,  x.  In  that  case,  the  sampled  interferogram  will  be  [15]; 

I'(x)  =  I(x)T(x)^(x  +  6)  (42) 

where  5  is  the  phase  error  in  'F,  and  T(x)  is  the  truncation  function  used  to  account  for  the 
fact  that  the  mirror  is  scanned  over  a  finite  path  difference,  instead  of  the  infinite  path 
difference  assumed  in  equation  (35). 

Performing  a  cosine  Fourier  transform  on  equation  (42)  yields; 

B'(a)  =  [B(a)  *  t(a)]  (43) 
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where  B(a)  is  the  true  spectrum,  and  ^(a)  is  a  phase  error  resulting  from  an  error  in  the 
zero  path  difference  position  in  the  sampling  function,  as  well  as  from  imperfect 
compensation  for  the  beamsplitter. 

The  consequence  of  the  phase  errors  include  distortion  in  the  spectral  lines,  as 
well  as  a  modification  of  the  baseline  of  the  spectrum.  Methods  of  elimination  of  the 
phase  error  are  discussed  below. 

4.  Elimination  of  Phase  Error.  If  the  phase  function  (j)(cr)  is  known,  the  sampled 
interferogram  can  be  corrected  by: 

l"(x)  =  I(x)  *  (44) 

In  the  Bomem  Fourier  transform  spectrometer,  the  phase  function  may  obtained  by 
recording  a  low-resolution  interferogram  on  both  sides  of  zero  path  difference.  The  phase 
function  is  calculated  from  [16]: 

(l)(o)  =  -  arctanj  (45) 

where  mi(CT)  and  mr(cy)  are  respectively  the  imaginary  and  real  components  of  the  Fourier 
transform  of  the  interferogram  as  a  function  of  o. 

5.  Limits  of  Resolution.  If  tilt  and  synchronization  errors  are  negligible,  an 
expression  for  the  inherent  resolution  of  the  spectrometer  may  be  derived  [15]. 
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The  number  of  spectral  elements  may  be  expressed  as: 


m 


(46) 


where  the  spectrum  ranges  from  ai  to  02 ,  and  5c  is  the  resolution  of  the  spectrometer. 

If  the  minimum  resolution  occurs  when  two  nearby  lines  are  resolved  for  a  path 
difference  L,  boundary  conditions  require: 


L 


—  and  L 
c 


m  +  1 
c  +  5c 


(47) 


Eliminating  L  from  the  above  equation  yields: 

8a  =  5/  (48) 

Expressing  m  in  terms  of  the  path  difference  gives  an  order  of  magnitude  estimate 
for  the  minimum  resolution: 


5c  «  (49) 

In  practical  situations,  the  minimum  achievable  resolution  for  the  instrument  is 
larger  than  this.  Thus  the  resolution  of  the  instrument  may  be  expressed  in  terms  of  the 
resolving  power,  R  =  c/Ac,  where  R  is  given  by  the  manufacturer,  and  c  is  an  average 
wavenumber  at  which  the  spectra  is  recorded.  R  will  depend  on  several  factors,  with 
mechanical  stability  and  solid  angle  subtended  by  source  being  dominant. 
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III.  Equipment  and  Software 


A.  Discharge  Tube 

The  discharge  tube  is  U-  shaped  (see  Figure  6),  and  mounted  upright  and  level, 
with  the  left  pipe  nearest  to  the  spectrometer.  The  hollow  electrodes  are  mounted 
vertically  in  the  sidearms.  The  sidearms  are  15  cm  in  length  (measured  from  electrode 
surface  to  the  near  wall  of  horizontal  section),  and  the  horizontal  section  of  the  tube  is 
approximately  9.5  cm  in  length,  with  an  sidearm  separation  of  about  5  cm.  The  inner 
diameter  of  the  tube  is  0.8  cm.  Optical  windows  are  mounted  on  the  ends  of  the 
horizontal  section. 


The  hollow  electrodes  and  mounting  rods  were  machined  from  Kovar.  The  inner 
diameter  of  the  hollow  electrodes  is  2.54  cm,  and  was  hollowed  to  a  depth  of  3.810  cm. 
The  wall  thickness  of  the  electrodes  was  0.312  cm.  All  electrode  surfaces  were  polished. 
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Four  tungsten  probes  were  inserted  into  the  tube  for  measurements  of  the  electrical 
potential  relative  to  ground  (see  Figure  7).  The  joint  of  the  probe  housing  with  the  wall 
was  funnel-shaped,  leading  to  a  constriction.  The  funnel  shape  avoids  sharp  edges  that 
would  affect  the  boundary  conditions  of  the  plasma,  while  the  constriction  effectively 
reduces  the  surface  area  of  the  wire  exposed  to  the  plasma,  eliminating  the  need  to  insulate 
the  wire.  The  tips  of  probes  in  the  horizontal  section  extended  1  mm  from  the  wall  of  the 
tube,  while  the  probes  in  the  side  arms  extended  to  the  center  of  the  tube. 


The  tube  was  filled  by  flowing  gas  through  the  left  pipe  of  the  discharge  tube.  For 
the  static  gas  fill  experiments,  the  right  pipe  was  sealed,  and  the  tube  was  pumped  out 
through  the  left  pipe.  For  the  flowing  gas  experiments,  the  gas  entered  the  tube  from  the 
left  pipe,  and  was  pumped  from  the  right  pipe.  Details  on  the  vacuum  system  are  included 
in  the  next  section. 
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B.  Vacuum  System  (Including  Flow  Svstem't 

Two  separate  vacuum  system  configurations  were  used  for  this  experiment.  The 
first  configuration  was  for  a  static  gas  fill,  while  the  second  configuration  was  for  the 
flowing  gas  experiment. 

I.  Static  Gas  Fill.  For  the  static  gas  fill  configuration,  a  custom-built  pump 
station  was  used  (see  Figure  8).  This  was  an  oil-free  system,  using  an  Alcatel  CFV  100 
turbopump,  an  Alcatel  CFV  10  drag  pump,  and  a  KNF  oil-fi'ee  diaphragm  pump. 
Butterfly  valves  placed  at  the  vacuum  pump  ports  could  be  used  to  isolate  the  vacuum 
pumps. 

The  connection  between  the  pumpstation  and  the  discharge  tube  was  via  a  stainless 
steel  bellows.  The  bellows  was  sealed  to  the  glass  tube  via  an  0-ring.  In  addition,  the 
discharge  tube  could  be  isolated  from  the  pumpstation  with  a  stopcock. 

The  gas  was  supplied  through  a  manifold  which  could  be  isolated  using  Varian 
viton-sealed  vacuum  valves,  and  a  needle  valve  at  the  gas  bottle.  Connections  were 
available  on  the  manifold  for  1  liter  glass  flasks,  as  well  as  lecture  bottles.  Matheson 
research  grade  helium  (99.9999%  pure)  and  Spectra  Gases  research  grade  nitrogen 
(99.9995%  pure)  were  used  for  this  phase  of  the  experiment. 

Pressure  was  measured  at  the  pumpstation,  using  a  10  torr  MKS  Baratron,  and  a 
100  torr  Vacuum  General  pressure  transducer.  In  addition,  a  Granville-Phillips  ion  gauge 
was  used  to  measure  pressures  <10"^  torr. 
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2,  Flowing  Gas  System.  The  flowing  gas  system  can  be  broken  down  into  three 
functions:  mass  flow  control,  pump  rate  control,  and  pumping. 

The  gas  mass  flow  rate  was  controlled  using  MKS  type  1259C  mass  flow 
controllers.  The  flow  rates  were  set  and  monitored  using  an  MKS  type  247  C  flow 
controller  readout.  For  this  portion  of  the  experiment,  Airco  Grade  6  (99.9999%  pure) 
helium  and  Matheson  UHP  (99.999%  pure)  nitrogen  were  used.  Pressure  was  measured 
upstream  from  the  discharge  tube  using  an  MKS  Baratron  pressure  transducer. 


An  Alcatel  2012A  mechanical  pump  was  used  for  the  flowing  gas  experiments. 
The  pump  rate  was  controlled  by  an  MKS  type  253A  butterfly-type  control  valve,  placed 
in  series  with  the  vacuum  pump.  The  butterfly-valve  was  controlled  by  an  MKS  type 
252A  Exhaust  Valve  Controller,  which  monitored  the  pressure  measured  by  the  Baratron, 
and  generated  an  error  signal  to  control  the  butterfly-valve.  The  mechanical  pump  could 
be  isolated  with  a  Varian  viton-sealed  vacuum  valve.  Overnight  pumpdown  was 


28 


accomplished  using  the  turbopump,  drag  pump,  and  diaphragm  pump  system  described 
above. 


C.  DC  Discharge  Power  Supply. 

The  discharge  was  powered  by  a  Glassman  High  Voltage  DC  power  supply  (model 
PS/EW10R60.0-11),  operating  in  the  current-controlled  mode.  The  current  was 
controlled  via  the  setting  of  a  potentiometer  on  the  power  supply  front  panel.  Current 
regulation  for  this  power  supply  is  rated  at  better  than  0.05%  [17].  The  accuracy  of  the 
current  display  is  rated  at  ~  1%. 


D.  Fourier-Transform  Spectrometer. 

A  Bomem  DA3.002  Fourier-transform  spectrometer  was  used  to  collect  the 
spectra.  The  essential  features  of  this  system  are;  the  input  optics,  a  Michelson 
interferometer,  a  detector,  and  a  vector  processor  (see  Figure  9). 

The  input  optics  consist  of  the  collection  optics,  an  adjustable  iris  aperture,  and  a 
collimator.  The  collection  optics  are  used  to  focus  an  image  of  the  source  onto  the  iris 
aperture.  The  collection  optics  consist  of  a  lens,  a  flat  steering  mirror,  and  an  off-axis 
paraboloid  mirror.  The  focal  length  of  the  lens  was  approximately  13  cm,  and  the  focal 
length  of  the  paraboloid  mirror  was  12.8  cm.  This  resulted  in  a  magnification  of  the  image 
of  the  tube  on  the  input  aperture  of  approximately  1:1. 
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The  iris  aperture  is  the  independent  factor  determining  the  resolving  power  of  this 
spectrometer,  as  it  limits  the  solid  angle  subtended  by  the  source.  The  aperture  was 
manually  set  by  closing  the  iris  around  a  1  nun  gauge.  Using  an  equation  supplied  by  the 
manufacturer,  the  resolving  power  is  limited  to  4.2  x  10^  for  an  aperture  diameter  of 
1  mm. 

The  Michelson  interferometer  was  used  to  generate  the  interferograms,  from  which 
the  spectra  could  be  obtained.  The  maximum  path  difference  for  this  model  was  250  cm. 
The  mirrors  were  dynamically  aligned,  using  electronic  tilt  correction,  resulting  in  less  than 
1  pradian  rms  tilt  error  over  the  full  travel  of  the  moviing  mirror.  The  velocity  accuracy  of 
the  moving  mirror  is  rated  at  ±  0.5%,  with  a  uniformity  of  ±  2%  or  +  0.002  cm/sec 
(whichever  is  greater). 

A  single  mode,  thermally  stabilized  He-Ne  laser,  operating  at  632.8  nm,  was 
coaxially  aligned  with  the  input  beam  of  the  interferometer.  The  interferogram  of  the  laser 
is  nearly  a  cosine  wave,  and  is  used  to  produce  sample  pulses  to  trigger  the  analog  to 
digital  conversion  of  the  interferogram.  In  normal  operation,  a  removable  mask  is  placed 
to  block  the  output  of  the  laser  to  the  detector.  This  mask  was  removed  for  this 
experiment  to  allow  light  propagating  along  the  axis  to  reach  the  detector.  The  light  from 
the  laser  was  blocked  from  the  detector  using  a  narrowband  filter. 

In  addition  to  the  He-Ne  laser,  a  broadband  white-light  source  is  coaxially  aligned 
with  the  input  beam.  The  purpose  of  the  white-light  source  is  to  produce  a  strong  signal 
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at  zero-path  difference  (determining  its  location),  which  then  allows  interferograms  from 
multiple  successive  scans  to  be  co-added. 

The  detector  was  an  R928  photomultiplier  (PMT),  manufactured  by  Hamamatsu 
[18],  It  has  a  nine-stage,  circular-cage  design.  The  input  window  material  is  UV  glass, 
and  the  photocathode  material  is  multialkali.  The  R928  has  a  spectral  range  of  285-900 
nm  with  a  peak  wavelength  at  400  nm. 

The  voltage  signal  from  the  detector  is  converted  from  analog  to  digital,  and 
stored  in  the  vector  processor.  The  vector  processor  co-adds  interferograms  from 
successive  scans  of  interferometer.  At  the  completion  of  the  scans,  the  vector  processor 
performs  a  fast-Fourier  transform  (FFT)  on  the  interferogram  to  obtain  the  spectrum.  The 
resulting  spectrum  and  interferogram  are  sent  over  a  local  network  to  a  VAX  computer. 


Figure  9.  Schematic  of  Fourier  Transform  Spectrometer 
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E.  Electric  Field  Measurement  System. 


As  described  in  section  A  of  this  chapter,  tungsten  probes  were  used  to  measure 
the  electric  potential  relative  to  ground  at  four  locations  in  the  positive  column.  Voltages 
were  recorded  by  an  automated  system  (see  Figure  10),  consisting  of  two  Fluke  model 
80K-40  HV  1000:1  voltage  reduction  probes.  The  resistance  of  the  Fluke  probes  was 
1000  MD,  (a  high  resistance  being  essential  to  minimize  current  drawn  from  the 
discharge).  The  probe  voltages  were  measured  with  two  HP3478A  multimeters.  The 
voltage  data  was  sampled  with  a  HP9836C  computer  which  acquired  the  data  over  an 
IEEE  interface.  The  computer  used  a  program  written  in  BASIC  to  compute  the  mean 
and  standard  deviation  of  the  voltage. 

The  distances  from  the  left  electrode  to  each  tungsten  probe  were  measured  prior 
to  the  experiment,  allowing  the  electric  field  to  be  determined  from  the  slope  of  the 
voltage  plotted  as  a  function  of  distance. 


Figure  10.  Electric  Field  Measurement  System 
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F.  Peak  Picking  Programs. 

Two  routines  were  used  and  compared  for  determining  peak  intensity  wavelengths 
for  the  spectra.  The  first  routine  is  BOMFIT,  developed  in-house  at  Wright  Laboratory. 
BOMFIT  uses  a  least-squares  approach  to  fit  Doppler  lineshapes  to  the  data,  using  line 
position  and  temperature  as  adjustable  parameters  (instrument  lineshape  also  factored  into 
the  fit).  The  routine  used  the  Levenberg-Marquardt  method  to  vary  the  parameters  until 
the  least-squares  error  between  the  calculated  lineshape  and  the  experimental  lineshape 
was  minimized.  Data  output  from  the  program  include  line  positions,  line  intensities,  and 
temperature. 

The  second  routine  used  for  peak  finding  was  the  differential  peak  finding  function 
included  with  the  Bomem  spectrometer  software.  This  function  numerically  differentiates 
the  spectrum.  Sharp  oscillations  in  the  differentiated  spectrum  correspond  to  line 
positions.  A  threshold  level  is  set  using  a  cursor,  and  peaks  are  assigned  to  wavenumbers 
where  the  oscillation  exceeds  the  threshold  level. 

The  performance  of  the  two  systems  was  comparable.  BOMFIT  showed  better 
performance  in  noisy  regions  of  the  spectrum.  The  least-squares  method  may  have  been 
more  immune  to  noise  than  the  differential  method.  In  addition,  BOMFIT  yielded  values 
for  the  ion  temperatures. 
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IV.  Procedure 


A,  Measurement  of  Tube  Dimensions 

Prior  to  mounting  the  discharge  tube,  the  tube  dimensions  were  measured  using  a 
metric  tape  measure.  The  dimensions  of  the  tube  that  were  measured  included  tube 
diameter,  length  of  horizontal  section,  side  arm  length,  side  arm  separation,  probe 
location,  and  probe  separation. 

For  each  dimension,  nine  measurements  were  taken,  and  a  mean  computed.  The 
error  in  the  measurement  was  taken  to  be  1mm  or  the  standard  deviation  (whichever  was 
greater). 

The  distance  between  the  tungsten  probes  and  the  left  electrode  was  defined  to  be 
the  length  of  the  line  running  along  the  axis  of  the  tube  from  the  electrode  lip  to  the  probe. 
This  approximation  will  tend  to  overestimate  the  probe  distance,  but  for  long  distances, 
the  error  is  small. 


B.  Calibration  of  FTS  System:  Phase  File 

As  mentioned  in  the  theory  section,  determination  of  the  phase  function  of  a 
spectrometer  allows  compensation  for  phase  errors  in  the  spectra.  The  phase  function  was 
obtained  by  recording  the  interferogram  of  a  broadband  lamp  on  both  sides  of  the  ZPD 
location  of  the  moving  mirror.  The  spectrometer  software  computes  the  sine  and  cosine 
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Fourier  transforms  from  the  interferogram,  and  produces  the  phase  function  using 
equation  (45). 


C.  Alignment  of  Discharge  Tube 


The  tube  was  mounted  on  an  X-Y  translation  stage  with  vertical  and  rotational 
adjustments.  Alignment  of  the  optical  system  was  done  using  a  three  separate  tests.  The 
first  test  was  to  align  the  image  on  the  front  and  rear  windows,  produced  by  backlighting 
the  spectrometer  with  a  white  light  source  (see  Figure  1 1).  The  image  sizes  and  locations 
on  the  front  and  back  windows  of  the  discharge  tube  were  obtained  using  a  white  index 
card.  The  tube  orientation  was  adjusted  such  that  the  white  light  image  was  centered  on 
both  the  front  and  back  windows,  with  the  front  and  back  images  being  of  equal  size. 


Figure  11.  Backlight  Alignment  Technique 
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The  second  alignment  test  was  to  turn  on  the  discharge,  and  observe  its  image  on 
the  spectrometer  aperture.  Proper  alignment  is  indicated  by  having  the  pupil  of  the 
aperture  in  the  center  of  the  image  of  the  discharge  tube. 

The  third  alignment  test  is  performed  using  the  spectrometer’s  co-axial  He-Ne 
laser  alignment  beam.  The  interferometer  reflects  a  portion  of  the  laser  beam  out  the 
spectrometer  aperture.  As  this  beam  is  used  for  aligning  the  spectrometer,  it  can  be  used 
to  check  the  alignment  of  the  tube.  This  is  not  a  precise  test,  as  the  He-Ne  laser  is  not 
precisely  co-axial  with  optical  axis,  but  instead  is  parallel  to  the  axis.  However,  it  is  usetul 
for  detecting  gross  misalignments  of  the  tube. 


D.  Acquisition  of  Spectra 

The  initial  step  in  the  acquisition  of  spectra  is  preparation  of  the  gas  mix.  For  the 
static  fill  case,  the  sequence  begins  by  pumping  the  discharge  tube  until  it  reaches  a  base 
pressure  <10'^  torr.  The  vacuum  pump  and  discharge  tube  are  isolated  from  the  rest  of 
the  vacuum  system  while  the  gases  are  premixed  at  the  pumpstation.  Gases  are  premixed 
using  the  law  of  partial  pressures  to  determine  the  ratio  of  N2  to  He.  After  premixing,  the 
valve  to  the  discharge  tube  is  opened,  and  the  gas  mix  enters  the  discharge  tube.  The  final 
pressure  is  read  with  the  Baratron  at  the  pumpstation.  The  ratio  of  the  pressure  in  the 
isolated  pumpstation  to  the  final  pressure  is  a  constant,  and  can  be  used  to  compute  the 
necessary  pressure  during  premix  to  obtain  a  desired  pressure  in  the  discharge  tube. 
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For  the  flowing  gas  case,  the  sequence  begins  with  setting  the  flow  setpoints  of  the 
flow  controllers  and  the  pressure  setpoint  of  the  exhaust  valve  controller.  Next,  the  flow 
controllers  are  filled  with  gas,  and  after  a  warmup  period,  the  flow  controllers  are  zeroed. 
Then  the  turbopump,  drag  pump,  and  diaphragm  pump  are  isolated  from  the  pumpstation. 
Gas  flow  is  initiated  at  this  point,  and  after  the  pressure  rises  above  1  torr,  the  valve  to  the 
mechanical  pump  is  opened. 

The  mass  flow  parameters  were  set  to  maintain  a  constant  residence  time  of  0.026 
seconds  for  the  pressure  range  over  which  this  experiment  was  performed.  The  limiting 
factors  were  the  resolution  of  the  20  seem  flow  controller,  and  the  pump  rate  through  the 
discharge  tube.  This  limited  our  experimental  conditions  to  pressures  above  0.5  torr  and 
0.08  seem  of  N2.  The  flow  rate  of  He  was  adjusted  accordingly  to  achieve  the  desired 
total  pressure  in  the  tube,  yet  maintain  the  constant  residence  time. 

Preparations  are  made  to  make  voltage  measurements  (see  section  F,  this  chapter), 
the  power  supply  is  coimected  to  the  anode  of  the  discharge  tube,  and  the  cathode  is 
connected  to  ground.  The  power  supply  is  turned  on,  initiating  the  discharge,  and  the 
current  level  is  set  using  the  potentiometer  on  the  front  panel.  The  voltage  on  the  PMT  is 
adjusted  such  that  the  analog  to  digital  converter  on  the  spectrometer  will  read  between 
15%  and  80%  of  the  scale. 

At  this  point,  the  menu  driven  Bomem  software  is  used  to  control  the  acquisition 
of  spectra.  The  operator  selects  the  desired  resolution  and  the  number  of  scans  taken  by 
the  interferometer.  In  this  experiment,  at  least  nine  scans  were  recorded  and  co-added  to 
generate  the  interferograms  for  every  reported  spectrum.  In  most  cases  sixteen  scans 
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were  recorded.  The  desired  resolution  was  set  at  0.06  cm  *  for  all  reported  spectra.  The 
time  necessary  to  obtain  a  spectra  under  these  conditions  was  approximately  40  minutes. 


E.  Determination  of  Ion  Temperatures 

The  ion  translational  and  rotational  temperatures  of  the  N2^  ion  were  determined 
from  the  spectra.  The  translational  temperature  was  obtained  from  the  BOMFIT 
software,  which  fit  Doppler  line  shapes  to  the  lines  in  the  observed  spectra,  with 
temperature  being  one  of  the  parameters. 

The  rotational  temperature  can  be  determined  from  [19:126]: 

lem  =  (J’  +  J"  +  1)  exp[-BJ'(J'  +  l)hc  /  kT]  (50) 


where  lem  is  the  emission  intensity,  Cem  is  depends  on  the  change  in  dipole  moment  and  the 
total  number  of  molecules  in  the  initial  vibrational  level,  B  is  the  fundamental  molecular 
rotational  constant  for  N2^,  Qr  is  the  sum  of  the  Boltzmann  factors  over  all  of  the 
rotational  states,  and  J’  and  J”  are  the  rotational  quantum  numbers  for  the  transition. 

In  practice,  these  terms  have  been  previously  calculated  [20]  giving  a  new  fitting 
equation  of: 


I 


em 


(51) 
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where  LT  and  ET  are  respectively  the  linear  and  exponential  terms  that  correspond  to  a 
particular  line  in  the  spectra,  and  Trot  is  the  rotational  temperature.  At  this  point,  the 
rotational  temperature  is  found  by  performing  a  least-squares  fit. 


F.  Voltage  Measurements 

At  the  beginning  of  a  laboratory  session,  the  HP9836C  computer  was  loaded  with 
the  voltage  data  acquisition  program,  and  given  an  initial  time/date  setting.  Prior  to 
recording  a  spectra,  tungsten  probes  1  and  4  (located  in  the  sidearms)  were  connected  to 
the  Fluke  HV  voltage  probes.  After  the  discharge  was  initiated,  the  voltage  was  sampled 
at  the  two  probe  locations.  When  the  spectrometer  finished  scanning,  voltages  at  probes  1 
and  4  were  sampled  again  and  checked  against  the  initial  reading  for  drift.  The  power  to 
the  discharge  was  disconnected,  and  the  tube  was  grounded  for  safety.  The  Fluke  voltage 
probes  were  then  connected  to  tungsten  probes  2  and  3  (in  the  horizontal  section).  The 
power  to  the  discharge  was  turned  on,  and  the  voltages  were  sampled  at  probes  2  and  4. 
For  each  probe  voltage  measurement,  the  voltage  was  sampled  256  times  over  period  of 
approximately  two  minutes. 
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A.  Data 


1 .  Static  Gas  Fill 


Electric  Field.  The  axial  electric  field  was  measured  using  probes,  as 
discussed  previously.  At  low  current  (1  mA),  the  measured  field  shows  poor  agreement 
with  Hong  and  Miller  (see  Figure  12),  especially  at  low  pressure,  and  diverges  from 
Saykally’s  measurement  at  6  torr.  The  trend  is  essentially  a  flat  line,  showing  no  pressure 


dependence  for  the  electric  field. 


Electric  Field  vs  Pressure  (  Static,  I  =  1  mA  ) 


•  Kee 
A  Miller 

O  Corrected  Miller 
XSaykally 


Figure  12.  E  vs  P,  Static  Fill,  1  =  1  mA 


The  electric  field  plot  changes  considerably  as  the  discharge  current  is  increased  to 
5  mA  (see  Figure  13).  At  pressures  <  3  torr,  the  measured  field  agrees  with  the  values 
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obtained  by  Hong  and  Miller  using  probes,  and  parallels  their  corrected  field  values. 
However,  at  pressures  >  3  torr,  the  measured  field  is  close  to  both  the  Hong  and  Miller 


corrected  field,  as  well  as  Radunsky  and  Saykally’s  value  obtained  from  velocity- 
modulation  laser  absorption  spectroscopy  [7], 


Electric  Field  vs  Pressure  (  Static,  I  =  5  mA ) 


•  Kee 
A  Miller 

O  Corrected  Miller 
X  Saykally 


Figure  13.  E  vs  P,  Static  Fill,  1  =  5  niA 


Similar  behavior  occurs  at  10  mA.  The  trend  is  to  follow  the  uncorrected  Hong 
and  Miller  data  at  low  pressures,  and  agree  with  their  corrected  data,  and  with  Saykally  at 
higher  pressures.  The  cross-over  point  for  this  behavior  appears  to  be  ~  1 .5  torr. 
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Electric  Field  vs  Pressure  (  Static,  1  =  10  mA  ) 
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Figure  14.  E  vs  P,  Static  Fill,  I  =  10  mA 


Doppler  Shift.  In  the  case  of  the  static  gas  fill,  the  measured  Doppler  shift 
differed  considerably  from  that  obtained  by  Hong  and  Miller  (see  Figures  15, 16  and  1 7). 
The  measured  shift  is  slightly  lower  than  the  Hong  and  Miller  data  at  low  pressures,  but 
diverges  widely  at  higher  pressures.  The  measured  peak  Doppler  shift  occurs  at  ~  1  torr, 
in  comparison  to  ~  2  torr  for  Hong  and  Miller. 
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Doppler  Shift  I  Doppler  Shift 


Doppler  shift  vs  Pressure  (  Static,  I  =  1  mA  ) 
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75.  Doppler  Shift  vs  Pressure,  Static  Fill,  7=7  mA 


Doppler  Shift  vs  Pressure  (  Static,  I  =  5  mA ) 
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Figure  16.  Doppler  Shift  vs  Pressure,  Static  Fill,  1  =  5  mA 


Temperature.  The  measured  N2^  rotational  and  translational  temperatures 
are  given  in  Figures  18,  19,  and  20.  The  temperatures  appear  to  be  independent  of 
pressure  and  electrode  orientation.  The  rotational  temperature  shows  a  moderate  level  of 
dispersion  at  a  discharge  current  of  1  mA,  with  the  dispersion  tightening  as  the  discharge 
current  is  increased.  Both  the  rotational  and  translational  temperatures  increase  with  the 
discharge  current  (see  Figure  21),  but  the  functional  dependence  of  the  increase  cannot  be 
determined  from  the  three  data  points  from  this  experiment,  and  the  data  point  reported  by 
Radunsky  and  Saykally  [21], 

One  noteable  feature  of  these  plots  is  the  ~  200  K  difference  in  the  translational 
and  the  rotational  temperatures,  implying  nonequilibrium  between  the  translational  and 
rotational  energies.  Similar  observations  were  reported  by  Radunsky  and  Saykally,  who 
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Figure  1 7.  Doppler  Shift  vs  Pressure,  Static  Fill,  I  -  10  mA 


Temperature.  The  measured  rotational  and  translational  temperatures 
are  given  in  Figures  18,  19,  and  20.  The  temperatures  appear  to  be  independent  of 
pressure  and  electrode  orientation.  The  rotational  temperature  shows  a  moderate  level  of 
dispersion  at  a  discharge  current  of  1  mA,  with  the  dispersion  tightening  as  the  discharge 
current  is  increased.  Both  the  rotational  and  translational  temperatures  increase  with  the 
discharge  current  (see  Figure  21),  but  the  functional  dependence  of  the  increase  cannot  be 
determined  from  the  three  data  points  from  this  experiment,  and  the  data  point  reported  by 
Radunsky  and  Saykally  [21]. 

One  notable  feature  of  these  plots  is  the  ~  200  K  difference  in  the  translational  and 
the  rotational  temperatures,  implying  nonequilibrium  between  the  translational  and 
rotational  energies.  Similar  observations  were  reported  by  Radunsky  and  Saykally,  who 
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reported  a  translational  temperature  of  1200  K  and  a  rotational  temperature  of  412  K  for  a 
discharge  current  of  120  mA  [21], 


Ion  Temperature  vs  Pressure  (  Static,  I  =  1  mA ) 
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. mean  rotational,  T  =  386  K 

—  -  —  -  mean  translational,  T  =  597  K 


Figure  18.  Ion  Temperature  vs  Pressure,  Static  Fill,  1=1  mA 


Ion  Temperature  vs  Pressure  ( Static,  I  =  5  mA ) 
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Figure  19.  Ion  Temperature  Pressure,  Static  Fill,  1  =  5  mA 
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Temperature  (Kelvin)  |  ^  I  Temperature  (Kelvin) 


Ion  Temperature  vs  Pressure  (  Static,  I  =  10  mA ) 
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20.  Ion  Temperature  vs  Pressure,  Static  Fill,  1=10  mA 
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2.  Flowing  Gas 

Electric  Field.  The  measured  electric  field  at  1mA  discharge  current  with 
a  flowing  gas  differs  considerably  from  the  case  of  the  static  gas  fill  (see  Figures  12  and 
22).  Comparing  the  two,  the  static  gas  fill  displayed  essentially  no  dependence  of  the  E 
field  due  to  pressure,  whereas  in  the  case  of  the  flowing  gas,  the  electric  field  is  pressure 
dependent.  It  is  observed  that  the  measured  field  closely  follows  the  uncorrected  Hong 
and  Miller  data  at  low  pressures,  and  parallels  their  corrected  data.  At  higher  pressure, 
the  measured  field  agrees  with  Radunsky  and  Saykally,  and  the  corrected  Hong  and  Miller 
data.  The  cross-over  point  for  the  behavior  is  ~  1 .5  torr. 


Electric  Field  vs  Pressure  (  Flowing  Gas,  I  =  1  mA ) 


Pressure  (torr) 


Figure  22.  E  vs  P,  Flowing  Gas,  1=1  mA 

At  discharge  currents  of  5  and  10  mA,  the  electric  field  displays  similar  behavior  to 
the  measurements  taken  at  1  mA  (see  Figures  23  and  24). 
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E  field  (volts/cm)  1^1  ^  (volts/cm) 


Electric  Field  vs  Pressure  (  Flowing  Gas,  I  =  5  mA  ) 
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23.  E  V5  P,  Flowing  Gas,  I  =  5  niA 


Figure  24.  E  vs  P,  Flowing  Gas,  I  ^  10  mA 
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Doppler  Shift.  In  the  case  of  the  flowing  gas,  the  measured  Doppler  shift 
showed  similar  behavior  to  the  case  of  the  static  gas  fill.  At  low  pressures,  the  size  of  the 
shift  was  lower  than  the  Hong  and  Miller  data.  At  pressures  between  ~  0.75  to  ~  5  torr, 
the  shift  values  diverge,  until  at  higher  pressures,  they  agree.  The  peak  value  of  the 
Doppler  shift  occurred  at  ~  0.75  torr  for  discharge  currents  of  5  and  10  mA,  in 
comparison  to  ~  2  torr  for  Hong  and  Miller.  The  highest  value  of  the  Doppler  shift  at  1 
mA  occurred  at  the  lowest  pressure,  0.5  torr. 


Figure  25.  Doppler  Shift  vs  Pressure,  Flowing  Gas,  I  =  J  mA 
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Doppler  Shift  vs  Pressure  (  Flowing  Gas,  I  =  5  mA ) 


7.00E-07 

6.00E-07 

5.00E-07 

C 

M  4.00E-07 

i:  3.00E-07 

o 

a 

2.00E-07 
1  .OOE-07 
0.CX)E+00 


0  12  3  4 

Pressure  (torr) 


Figure  26.  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1  =  5  /nA 
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Doppler  Shift  vs  Pressure  (  Flowing  Gas,  1  =  10  mA ) 
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Figure  27.  Doppler  Shift  vs  Pressure,  Flowing  Gas,  I  =  10  mA 

Temperature.  The  rotational  and  translational  temperatures  are  given  in 
Figures  28,  29,  and  30  for  the  case  of  flowing  gas.  Significant  differences  are  observed  in 


50 


comparison  to  the  static  gas  fill.  The  mean  rotational  temperatures  are  lower  for  the 
flowing  gas  than  for  the  static  gas  fill.  This  is  to  be  expected,  as  residence  time  of  the 
flowing  gas  is  short  (  ~  0.026  sec),  and  thus  the  flowing  gas  can  be  expected  to  stay  at 
room  temperature.  In  contrast,  the  mean  translational  temperatures  are  slightly  higher  for 
the  flowing  gas  than  for  the  static  gas  fill. 

The  rotational  temperatures  show  a  wide  dispersal  at  a  discharge  current  of  1  mA, 
with  some  of  the  data  points  being  below  room  temperature.  This  is  attributed  to  low 
signal-to-noise  of  the  spectra.  This  dispersal  is  sharply  reduced  as  the  discharge  current  is 
increased.  Similar  behavior  is  observed  for  the  translational  temperature. 


Ion  Temperature  vs  Pressure  (  Flowing  Gas,  I  =  1  mA  ) 
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Figure  28.  Ion  Temperature  v.v  Pressure,  Flowing  Gas,  1^1  mA 
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Temperature  (Kelvin)  |  ^  I  Temperature  (Kelvin) 


Ion  Temperature  vs  Pressure  (  Flowing  Gas,  I  =  5  mA  ) 


29.  Ion  Temperature  vs  Pressure,  Flowing  Gas,  1  =  5  mA 


Figure  30.  Ion  Temperature  vs  Pressure,  Flowing  Gas,  I  =  10  mA 
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Temperature  (Kelvin) 


Ion  Temperature  vs  Current  (Flowing  Gas) 


Figure  31.  Ion  Temperature  vs  Current,  Flowing  Gas 
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B.  Analysis 

1.  Corrections  to  Electric  Field  in  Low  Pressure  Limit.  The  electric  fields 
resulting  from  probe  measurements  resulted  in  poor  agreement  with  previously  published 
data  at  low  pressures,  but  good  agreement  at  high  pressures.  To  correct  the  low  pressure 
data,  the  electric  field  was  calculated  from  the  measured  Doppler  shift  using  equation  (5). 
As  the  collision  rate  increases  with  pressure,  equation  (5)  is  no  longer  a  good 
approximation,  and  as  the  following  figures  show,  the  calculated  electric  field  decreases 
with  pressure  after  reaching  a  peak  at  ~  1  torr. 


The  adjusted  electric  field  will  then  consist  of:  1)  at  low  pressure  (  P  <  1  torr), 
values  calculated  from  the  measured  Doppler  shift;  2)  at  high  pressure  (3-6  torr),  probe 
measurements;  and  3)  linearly  interpolated  values  for  the  mid-pressure  range  (1-3  torr). 
The  adjusted  electric  fields  are  given  in  Figures  32-37. 


Figure  32.  Adjusted  Electric  Field  vs  Pressure,  Static  Fill,  1=1  niA 
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E  field  (volts/cm)  ^1  ^  (volts/cm) 


Electric  Field  vs  Pressure  (  Static,  I  =  5  mA  ) 


55.  Adjusted  Electric  Field  vs  Pressure,  Static  Fill,  1  =  5  mA 
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Figure  34.  Adjusted  Electiic  Field  vs  Pressure,  Static  Fill,  J  =  10  mA 
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Electric  Field  (volts/cm)  ^  |  Electric  Field  (volts/cm) 


Electric  Field  vs  Pressure  (  Flowing  Gas,  I  =  1  mA ) 


35.  Adjusted  Electric  Field  vs  Pressure,  Flowing  Gas,  1=1  mA 


Electric  Field  vs  Pressure  (  Flowing  Gas,  I  =  5  mA ) 
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Figure  36.  Adjusted  Electric  Field  vs  Pressure,  Flowing  Gas,  1  =  5  tnA 
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Electric  Field  vs  Pressure  (  Flowing  Gas,  1  =  10  mA  ) 
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Figure  37.  Adjusted  Electric  Field  V5  Pressure,  Flowing  Gas,  I  -  10  niA 


2,  E-scaled  Doppler  Shift  vs  Pressure.  The  E-scaled  Doppler  shifts  were 
obtained  by  dividing  the  measured  Doppler  shift  by  the  adjusted  electric  field.  The  results 
are  displayed  in  Figures  38-43.  In  general,  there  is  poor  agreement  between  the  measured 
values  and  the  previously  published  values  of  Hong  and  Miller.  The  E-scaled  Doppler 
shift  tends  to  follow  the  same  trend  as  the  Hong  and  Miller  data,  but  the  expected  drop  in 
the  shift  due  to  collisions  takes  place  at  lower  pressures  than  in  their  data. 
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E-scaled  Doppler  shift  (cm/volt)  ^  I  E-scaled  Doppler  shift  (cm/volt) 


E-scaled  Doppler  shift  vs  Pressure  (  Static,  I  =  1  mA ) 
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38.  E-scaled  Doppler  Shift  V5  Pressure,  Static  Fill,  1  =  1  mA 
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E-scaled  Doppler  shift  (cm/volt)  |  ^  j  E-Scaled  Doppler  Shift  (cm/volt) 


E-Scaled  Doppler  Shift  vs  Pressure  (  Flowing  Gas,  I  =  5  mA ) 


42.  E-scaled  Doppler  Shift  vs  Pressure,  Flowing  Gas,  1  =  5  mA 
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3.  Disagreement  with  Hong-Miller  Experiment.  For  all  experimental  conditions, 
the  observed  Doppler  shift  disagreed  with  the  Hong-Miller  data.  This  experiment 
operated  under  similar  conditions  to  that  of  Hong  and  Miller  with  the  exception  that  the 
discharge  currents  were  lower  (1,5,  and  10  mA  vs  25  mA). 

A  possible  cause  for  the  disagreement  can  be  found  by  examining  the  four  likely 
processes  by  which  the  N2^(B)  ion  are  formed  (see  Figure  44).  The  first  likely  process  is  a 
one-step  process,  where  N2  (B)  is  created  by  collision  of  a  ~  19  eV  electron  with  N2.  The 
second  likely  process  is  a  two-step  process  where  first  N2  (X)  is  created  by  collision  of  a 
~  16  eV  electron  with  N2,  then  N2^(B)  from  the  collision  of  the  N2  (X)  ion  with  a  ~  3  eV 
electron.  The  third  and  fourth  processes  involve  collision  of  a  He  metastable  with  N2. 
The  metastable  ionization  can  leave  the  N2^  ion  in  either  the  B  or  X  state.  If  the  ion  is  left 
in  the  B  state,  it  rapidly  decays  into  the  X  state  via  radiative  transition  (  t  =  66ns).  Then 
N2^(B)  can  be  formed  by  collision  with  a  ~  3  eV  electron. 


Figure  44.  N2^  (B)  Formation  Processes 
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The  method  of  N2^(B)  formation  affects  the  observed  Doppler  shift.  The  theory 
developed  in  Chapter  2  assumes  a  one-step  process.  If  a  two-step  process  is  involved, 
then  the  observed  Doppler  shift  would  increase,  because  the  ion  accumulates  kinetic 
energy  from  the  applied  electric  field  while  in  the  X  state. 

Examination  of  the  electron  energy  distribution  fimction  for  the  two  experimental 
conditions  can  be  used  to  qualitatively  argue  that  the  Hong-Miller  experiment  may  be 
operating  under  conditions  that  include  a  greater  fraction  of  two-step  processes  than  in  the 
author’s  experiment.  In  Figure  45,  the  electron  energy  distribution  functions  are  plotted 
for  discharge  currents  of  10  and  25  mA.  The  plot  suggests  that  the  one-step  process  may 
compose  a  higher  fraction  of  the  total  rate  of  N2^(B)  formation  for  the  10  mA  case,  than 
for  the  25  mA  case.  However,  a  detailed  calculation,  beyond  the  scope  of  this  thesis,  is 
required  to  definitively  determine  the  relative  contributions  of  one-  and  two-step  processes 
to  the  total  N2^(B)  formation  rate. 


Electron  Energy  Distribution  Function  (  P  =  2  torr ) 


- 1  =  10  mA  (Kee) 

<0-  I  =  25  mA  (Miller) 


Figure  45.  Electron  Energy  Distribution  Function  [22] 
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For  the  case  of  the  metastable  processes,  the  electron  energy  distribution  function 
suggests  that  at  ~  3  eV,  the  higher  current  has  a  larger  distribution  of  electrons  that  will 
collisionally  excite  N2^(X)  to  N2^(B).  Again,  the  residence  time  in  the  X  state  will  show  as 
a  bias  towards  larger  Doppler  shifts. 

4.  Plot  of  Surface.  The  problem  of  fitting  an  equation  to  a  set  of  data  by 
choosing  optimum  fitting  parameters  has  been  addressed  by  Bevington  [24:96-167].  The 
method  used  in  analyzing  data  for  this  experiment  is  the  least-squares  method.  Consider 
a  set  of  j  data  points  (y{xj),  xj),  where  xj  is  the  independent  variable,  and  y(xj)  is  the 
dependent  variable.  A  fitting  equation  F(xj,  a,  ft)  is  proposed  to  graph  against  the  data, 
where  a  and  ft  are  adjustable  parameters,  chosen  to  ‘optimize’  the  fit  to  the  experimental 
data. 

The  figure  of  merit  for  the  fit  is  chi-squared,  which  is  defined  as: 

x"  =  E  (f(Xj,  a,  P)  -  y(X;))’  4  (52) 

j 

where  cj  is  the  experimental  error  associated  with  the  jth  data  point.  Optimization  of  the 
fit  is  performed  by  choosing  values  of  a  and  ft  that  minimize  the  value  of 

To  determine  Qv  and  Qq,  the  least-squares  method  is  used  to  fit  equation  (25)  to 
the  experimental  E-scaled  Doppler  shifts.  Initially,  a  MathCAD  5.0  spreadsheet  was 
written  to  perform  the  fit,  using  the  Levenberg-Marquardt  algorithm  contained  in  the 
‘solve  block’  routine  [25:499].  This  proved  to  be  unsatisfactory,  as  the  value  obtained 
from  the  fit  proved  to  be  highly  sensitive  to  initial  conditions. 
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Although  it  is  computationally  expensive,  the  grid-search  method  was  finally  used, 
due  to  its  robustness,  to  search  for  the  minima  of  x^-  In  this  method,  is  computed  at 
trial  values  of  the  fitting  parameters,  and  the  minima  is  found  by  inspection.  A  MathCAD 
spreadsheet  was  written  to  plot  x^  as  function  of  Qv  and  Qq,  and  to  display  the  values  in  a 
table.  The  plots  are  given  in  Figures  47-52,  and  x^  values  are  presented  in  Tables  1-6. 

The  results  are  disappointing,  in  that  the  plots  show  x^  being  a  very  weak  function 
of  Qv,  although  being  strongly  dependent  on  Qq.  No  minima  in  x^  is  observed  for  positive 
values  of  Qv,  which  eliminates  the  possibility  of  determining  the  ion  mobility  from  this 
experiment. 

A  rough  estimate  of  the  quenching  cross-section  Qq,  can  be  obtained  from  the  chi- 
square  plots.  Using  the  established  value  of  the  mobility  (Ko  =  19  cm^/volt/sec),  an 
estimate  for  the  collision-cross  section  can  be  obtained  from  the  Langevin  equation 
[26:162]; 


K„  =  0.815  X 


ok  +  M 


gas 


M;„ 


(53) 


where 


k 


1 

^gasQv 


(54) 


using  the  standard  temperature  (273  K)  and  pressure  (760  torr).  Assuming  an  ideal  gas 
results  in  a  collision  cross-section  of  ~  34.2  angstrom^.  (This  is  in  disagreement  with  the 
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Hong-Miller  value  of  17.8  angstrom^.  This  may  be  the  result  of  their  expression  for  the 
collision  frequency  being  in  error,  reference  equations  28  and  29). 

Inspecting  the  plots  shows  that  for  the  static  gas  fill  experiments,  a  Qq  value  on  the 
order  of  15-20  angstrom^  is  anticipated,  white  the  flowing  gas  experiment  results  in  a  Qq 
value  of  20-25  angstrom^.  In  comparison,  Hong  and  Miller  reported  a  value  of  6.3 
angstrom^.  The  higher  quenching  cross-section  for  the  flowing  gas  rate  is  attributed  to  a 
higher  ion-gas  collision  velocity  due  to  the  flowing  of  the  gas  through  the  tube.  The 
difference  between  our  measured  values  and  Hong  and  Miller’s  is  attributed  to  the  bias  in 
the  observed  Doppler  shift  which  results  from  a  significant  fraction  of  two-step  N2^(B) 
production. 

The  larger  quenching  rate  for  the  flowing  gas  in  comparison  to  the  static  gas  fill 
corresponds  to  a  reduced  measured  E-scaled  Doppler  shifts  in  the  mid-pressure  regions,  as 
well  as  a  shift  in  the  peak  Doppler  shifts  to  lower  pressures.  This  also  holds  for 
comparisons  to  the  Hong  and  Miller  data 
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Table  1.  ^  values,  Static  Fill,  1  =  1  mA 


a 

angstrom^ 

a 

angstrom^ 

11 

12 

13 

14 

15 

16 

1 

32.969 

24.936 

20.899 

20.259 

22.559 

27.435 

34.578 

2 

32.617 

24.847 

21.026 

20.564 

23.013 

28.014 

35.261 

3 

32.289 

24.779 

21.170 

20.884 

13A19 

28.603 

35.953 

4 

31.986 

24.732 

21.331 

21.217 

23.957 

29.202 

36.653 

5 

31.706 

24.703 

21.508 

21.564 

24.447 

29.81 

37.361 

6 

31.451 

24.695 

21.701 

21.925 

24.947 

30.428 

38.077 

7 

31.219 

24.706 

21.911 

22.299 

25.459 

31.055 

38.801 
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Figure  47.  ^  Plot,  Static  Fill,  I  =  5  mA 


Table  2.  ^  values.  Static  Fill,  I  =  5  mA 


0V 

angstrom^ 

Q. 

angstrom^ 

14 

15 

16 

17 

18 

19 

20 

1 

66.299 

51.093 

41.792 

37.662 

38.031 

42.292 

49.903 

2 

65.416 

50.604 

41.629 

37.769 

38.361 

42.807 

50.570 

3 

64.562 

50.137 

41.484 

37.892 

38.705 

43.333 

51.247 

4 

63.734 

49.693 

41.358 

38.029 

39.061 

43.870 

51.932 

5 

62.934 

49.271 

41.250 

38.182 

39.430 

44.417 

52.626 

6 

62.160 

48.871 

41.161 

38.151 

39.812 

44.975 

53.329 

7 

61.413 

48.493 

41.09 

38.534 

40.206 

45.543 

54.04 
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Table  3.  values,  Static  Fill,  1  =  10  mA 


Ov 

anf^strom^ 

a 

angstrom^ 

17 

18 

19 

20 

21 

1 

198.811 

191.984 

190.631 

193.901 

201.06 

2 

198.639 

192.103 

190.988 

193.450 

201.766 

3 

198.484 

192.237 

191.357 

195.010 

202.48 

4 

198.346 

192.385 

191.738 

195.579 

203.202 

5 

198.226 

192.548 

192.130 

196.158 

203.931 

6 

198.124 

192.725 

192.534 

196.746 

204.669 

7 

198.038 

192.915 

192.949 

197.344 

205.414 
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Table  4.  values,  Flowing  Gas,  I  =  I  mA 


Q. 

angstrom^ 

22 

23 

24 

25 

26 

1 

7.559 

7.193 

7.619 

8.756 

10.532 

2 

7.554 

7.228 

7.697 

8.872 

10.682 

3 

7.533 

7.266 

7.778 

8.991 

10.833 

4 

7.525 

7.307 

7.862 

9.112 

10.986 

5 

7.521 

7.351 

7.948 

9.235 

11.142 

6 

7.521 

7.398 

8.038 

9.361 

11.299 

7 

7.523 

7.449 

8.129 

9.489 

11.459 
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so 


Figure  50.  Plot,  Flowing  Gas,  1  =  5  mA 


Table  5.  values.  Flowing  Gas,  1  =  5  mA 


Qr 

angstrom^ 

Qo 

angstrom^ 

19 

20 

21 

22 

23 

1 

Mni>9 

43.340 

42.145 

43.667 

47.493 

2 

47.498 

43.282 

42.238 

43.885 

47.814 

3 

Ai.iei 

43.233 

42.339 

48.140 

4 

47.047 

43.192 

42.447 

48.472 

5 

46.838 

43.161 

42.563 

44.577 

48.808 

6 

46.640 

43.139 

42.686 

mSmm 

49.150 

7 

46.452 

43.126 

42.816 

49.497 
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Table  6.  values,  Flowing  Gas,  1  =  10  mA 


0v 

angstrom^ 

Q, 

angstrom^ 

22 

23 

24 

25 

26 

1 

29.511 

29.082 

29.385 

30.322 

2 

29.506 

29.177 

29.453 

30.420 

3 

30.685 

29.504 

29.154 

29.524 

30.519 

4 

30.639 

29.504 

29.193 

29.597 

30.620 

5 

30.596 

29.506 

29.234 

29.671 

30.722 

6 

30.555 

29.511 

29.278 

29.747 

30.826 

7 

30.518 

29.518 

29.323 

29.825 

30.932 

Comparison  to  Hong-Miller  results.  The  failure  to  obtain  a  physically 
realistic  collision  cross-section  via  this  experiment  raises  the  question  of  why  this 
experiment  failed  where  Hong  and  Miller  succeeded.  The  explanation  may  lie  in  a 
statistical  analysis  of  the  Hong-Miller  data.  For  this  purpose,  their  data  was  extracted 
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from  graphs  by  measuring  with  a  metric  ruler.  The  resolution  of  the  ruler  was  well  within 
the  8  mm  wide  error  bars  on  the  graph. 

The  values  of  the  fit  to  the  Hong-Miller  data  were  computed  and  are  presented 
below  in  Figure  53  and  Table  7.  The  plot  on  their  data  shows  a  weak  dependence  on 
Qv  and  a  strong  dependence  on  Qq.  Inspection  of  the  values  in  Table  7  yield  Qv  «  5 
angstrom^  in  comparison  to  their  published  value  of  17.8  ±  1.8  angstrom^,  and  Qq  »  7 
angstrom^  in  comparison  to  their  published  value  of  6.3  ±  0.7  angstrom^. 

This  large  disagreement  between  the  Qv  given  by  Hong  and  Miller,  and  the  values 
is  likely  due  to  the  flat  nature  of  the  y^  surface  along  the  Qv  axis.  The  Hong-Miller  data 
was  not  given  in  text  form,  nor  was  their  curve-fitting  method  specified.  Slight  deviations 
in  their  data,  and  the  values  extracted  from  the  graph  may  be  responsible  for  shifting  the 
minima  along  the  very  flat  y^  surface.  The  choice  of  curve-fitting  method  could  also 
impact  the  final  value  of  the  fitting  parameters.  For  example,  the  widely  used  Levenberg- 
Marquardt  method  performs  poorly  when  the  y^  surface  is  flat,  and  could  converge  to  an 
incorrect  result  [27:525]. 

The  apparent  sensitivity  of  the  location  of  the  minimum  to  either  a  small  change  in 
initial  conditions  (experimental  data),  or  a  change  in  the  minima  search  algorithm 
implies  that  the  1.8  angstrom^  uncertainty  in  Qv  given  by  Hong  and  Miller  is 
underestimated.  This  insensitivity  is  illustrated  in  Figure  54  for  the  case  where  equation 
(28)  is  used  to  define  the  collision  frequency. 
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Figure  52.  Plot,  Hong-Miller  Experiment,  I  =  25  mA 


Table  7.  ^  values,  Hong-Miller  Experiment,  I  =  25  mA 


0v 

angstrom^ 

Q, 

angstrom^ 

5 

6 

7 

8 

9 

1 

11.962 

4.543 

1.263 

1.152 

3.395 

2 

11.315 

4.252 

1.212 

1.263 

3.614 

3 

10.699 

3.983 

1.176 

1.384 

3.84 

4 

10.115 

3.735 

1.154 

1.514 

4.074 

5 

9.562 

3.506 

1.145 

1.654 

4.314 

6 

9.038 

3.298 

1.150 

1.803 

4.561 

7 

8.543 

3.109 

1.168 

1.962 

4.814 
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E-scaled  Doppler  Shift 


Figure  53.  Insensitivity  of  Fitting  Equation  to  Choice  of  Qv  for  Hong-Miller  Data. 


74 


VI.  Conclusions  and  Recommendations 


A.  Summary  of  Results 

This  experiment  failed  to  determine  the  mobility  of  the  N2^  ion  in  helium.  This  was 
due  to  the  inability  to  reduce  the  data  using  the  model  of  Hong  and  Miller.  A  surface 
was  calculated  from  a  fit  of  the  experimental  data  to  a  theoretical  expression  for  the  E- 
scaled  Doppler  shift.  This  surface  was  steep  in  the  direction  of  the  quenching  cross- 
section,  but  flat  in  the  direction  of  the  collision  cross-section.  Whether  the  fitting  equation 
is  inherently  insensitive  to  a  fit  to  the  collision  cross-section  has  not  been  proven.  It  is 
possible  that  the  flat  surface  is  the  result  of  the  model  not  being  applicable  to  the  data. 

Estimates  of  the  quenching  cross-section  were  obtained  from  the  plots, 
assuming  a  collision  cross-section  of  34.2  angstrom^  derived  from  the  Langevin  expression 
for  the  mobility,  and  previously  published  values  of  the  N2^  mobility  in  helium.  For  the 
static  gas  fill,  Qq  «  15-20  angstrom^,  while  for  the  flowing  gas,  Q,  »  20-25  angstrom^. 
This  is  higher  than  the  value  reported  by  Hong  and  Miller  (Qq  =  6.3  ±  0.7  angstrom^). 


B.  Recommendations  for  Further  Research 

A  statistical  analysis  of  the  data  from  this  experiment  and  the  Hong-Miller 
experiment  show  that  measurements  of  the  Doppler-shift  resulting  from  the  emission 
spectra  of  the  B  X  transition  of  N2^  yield  a  highly  uncertain  value  for  the  collision 
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cross-section.  Further  work  is  needed  on  the  model  presented  by  Hong  and  Miller  to 
determine  under  what  conditions  it  may  be  applied  towards  experimental  data 

This  leaves  open  the  question  of  switching  to  the  A  ->  X  transition  of  N2  . 
Whether  this  would  be  fruitful  depends  on  the  sensitivity  of  the  surface  to  the  fitting 
parameters.  Hong  and  Miller  have  published  cross-section  resulting  from  this  transition. 
A  statistical  analysis  of  their  results,  as  performed  here  for  the  B  — >  X  transition,  should  be 
accomplished  prior  to  performing  the  experiment.  In  circumstances  where  no  prior  data  is 
available,  a  rough  estimate  could  be  made  of  the  ion  mobility,  and  thus  the  collision-cross 
section.  Values  could  then  be  generated  for  the  E-scaled  Doppler  shift  (a  random  number 
generator  could  be  used  to  produce  statistical  error)  for  a  range  of  pressures. 
Examination  of  the  resulting  surface  would  then  reveal  whether  the  fitting  equation  is 
sensitive  to  the  collision  cross-section. 

Plots  of  the  electron  energy  distribution  function  for  experimental  parameters 
should  be  used  for  choosing  a  discharge  current  that  minimizes  the  Doppler  shift  bias 
associated  with  two-step  excited  ion  production  processes,  yet  allows  sufficient  signal  for 
recording  the  spectra.  Temperature  data  may  prove  useful  for  further  analysis  of  this 
diagnostic  technique.  For  this  reason,  a  thermocouple  should  be  used  to  measure 
discharge  tube  wall  temperatures.  Another  alternative  is  to  use  a  water  cooling  jacket 
around  the  discharge  tube  to  fix  the  wall  at  a  desired  temperature. 
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Appendix  A:  Example  Search  for  Minima  of 

Units: 

angstrom  :=10’^*’  m 

amu  :  =  1 .6605402 1 0'^’-kg  Atomic  mass  unit 


Constants: 

c  := 299792458—  Velocity  of  light 

sec 

k  D  :=  1.380658 10'^^-^—^  Boltzmann’s  constant 
®  K 

D  g  := 900  cm^-^^  reduced  ambipolar  diffusion  constant 
sec 


q  :=  1.60217733 10’*® 

coul  elementary  charge 

_9 

T  :=6610  sec 

Lifetime  of  B  State  of  N  2+  ion 

M  jjg  :=4  amu 

Mass  of  helium  atom 

M  :=28  amu 

Mass  of  N  2  molecule/ion 

M  He'^  N2 

H  := - 

M  He  +  M  n2 

|i  =  3.5 ‘amu 

reduced  collision  mass 
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Data: 

Ra  :=0.4  cm  Discharge  tube  radius 

T  :=597  K  Ion  translational  temperature 


Gas  Pressure 


E-scaled  doppler  shift 


7.58  lOT* 
7.58  lO'* 
7.5410'* 
4.3610'* 
1.8510'* 
1.8610'* 


cm 

volt 


.Ol  torr  Partial  pressure  of  nitrogen 


error  in  E-scaled 


3ler  shift 


®  Shift 


5.0110'^ 
2.45 10'^’ 
2.6510'® 
5.42 10'® 
3.410'® 
4.2510'® 


cm 

volt 


Helium  number  density 


relative  collision  velocity 
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M 


r(Qv,k) 


He 


M 


N2+-^He 


-  Qv  N  V 


Pa  T  /2.40^^ 
PHe'273KlRa  j 

k 


}^Qq,k) 


Nj-He  collision  frequency 


Wall  collision  rate 


Mean  free  path  for  quenchini 


Theoretical  expression  for  E-scaled  doppler  shift 


y(Qv,Qq,k) 


q 

1  -  exp 

1  »  ' 

1  X(Qq,k) 
i  V-t 

r((^',k)-)i(Qq,k)\ 
V  / 

MN2r(Qv,k)c 

i  +  r(Qv,k)-T 

1  -  expj 

1  ^Qq,k)^ 
1  V-t  1 

vel  :=1..50 
quench  :=1..20 

Qv^^i  := vel- angstrom^  Trial  values  for  momentum-transfer  cross-section 


uench 


:  =  quench  •  angstrom 


2 


Trial  values  for  quenching  cross-sections 


Chi-square  calculation 


Z  squared  ,  , 

^  vel,  quench 
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Qq  (angstrom^) 


^  squared 


Qy  (angstrorrf ) 


Chi-SQuare  values 


Qq  (angstronf) 


Z  squared 


8 

9 

10 

11 

12 

13 

14 

1 

64.633 

45.79 

33.018 

25.098 

21.191 

20.691 

23.13$ 

2 

63.524 

45.115 

32.675 

25.02  ' 

21.328 

21.006 

23.602 

3 

62.455 

44.47 

32.357 

24.962 

21.482 

21.3351 

24.078 

4 

61.423 

43.855 

32.063 

24.924 

21.652 

21.67^ 

24.565 

5 

60.428 

43.27 

31.793 

24.905 

21.839 

22.035 

25.063 

6 

59.469 

42.713 

31.547 

24.906 

22.042 

22.40^ 

25.572 

7 

58.547 

42.185 

31.324 

24.926 

22.26  ; 

22.788 

26.093 

Qy  (angstrom?) 


Appendix  B:  Residence  Time  in  Discharge  Tube 


units: 

3 

seem  :=  12.66667lorr-^5^ 
see 


Data: 


Total  Pressure 


Helium  Flow  rate 


Nitrogen  Flow  rate 


P  := 


0.08' 

r-  - 

4 

0.5 

0.08 

5.994 

0.75 

0.08 

8.019 

1 

0.08 

1.5 

2 

12.069 

•torr  ^  He 

16.119 

•seem 

0.08 

0.08 

24.21 

3 

0.08 

32.32 

5 

0.08 

40.42 

5 

0.08 

48.52 

6 

■seem 


r  :=  0.4  cm  discharge  tube  radius 


X  :  =  51.667  mm  length  of  horizontal  section  of  discharge  tube 


Expressions: 

A  .  2 

A  :=7r-r 

A  =  0.503  ‘cm^  cross-sectional  area  of  discharge  tube 


tot 


N2  +  ^  He 
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■seem 


Total  mass  flow  rate 


tot 


4.08 

6.074 

8.099 

12.149 

16.199 

24.29 

32.4 

40.5 

48.6 


The  flow  velocity  is  given  by: 


PA 


V  = 


2.056 

2.041 

2.041 

2.041 

2.041 

2.04 

1.633 

2.041 

2.041 


ttrsee 


''mean  :=mean(v) 


''  mean 


199.729— 

see 


Then  the  residence  time  in  the  tube  is; 


X 

X  := - 

'mean 


"t  =0.026  see 
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Units 

amu  :=  1.660540210'^’ kg 


atomic  mass  unit 


angstrom  :=  10 


Angstrom 


Constants 

c  := 29979245^^51- 
sec 

q  :=  1.6021773310' ^®coul 
M  jjg  ;=  4.0026amu 
Mn2  •=2  14.0067amu 

kn  1.38065810' 

®  K 

Tg  :=6610  ^  sec 


Velocity  of  light 
Electron  charge 
Mass  of  Helium  atom 
Mass  of  N2  molecule 

Boltzman  constant 
Lifetime  of  lv^+  B  state 


T  :=600K 

2 

Qy  :=  34  angstrom 


''BOXO 


c 

39 14  angstrom 


Gas  temperature 

N2*-He  collision  cross-section 

frequency  of  &  Sigma  (v=0)  to  )^sigma  (v=0)  transition 


Frequency 


Wavenumber 


VcoU(P) 


kg-T  [  7t  ^Mjq2  Muej 


collision  frequency  as  a 
function  of  pressure 
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Pressure  broadening  in  frequency. 


Av 


pressure 


(P) 


1 

2  k 


rB  / 


pressure 


Av  fP) 

pressure  ' 


Pressure  broadening  in  wavenumber. 


Resolution  of  spectrometer 


Rj  ;=4.3  I0-^ 


Resolving  power  of  Bomem  spectrometer 


Aainst(o,R) 

Bomem^®^  -  Ao  ^a,R  jj 


Comparison  of  linewidths  for  B-X  fO.1)  transition  as  a  function  of  Pressure. 

i  :=  10,20..  760  iterator 

P.  :=i  torr  Pressure  variable 

_'’B0X0 

•^BX' - 

c 

.555*10'’  *0111'’  wavenumber  of  transition 


press  .  '  pressure  (^i) 

Pressure  broadening 

dopp .  •=  ^  doppler  BX"^) 

Doppler  broadening 

instr.  Boinem(°Bx) 

Instrument  broadening 

P. 

Ptorr.  :=— 

'  torr 
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O.I 


Linewidth  vs  Pressure 


’ '  ■  Pressure 
Doppler 
Instrument 


As  is  seen,  at  subatmospheric  pressures,  the  Doppler  linewidth  dominates.  The  instrument 
broadening  is  comparable  to  the  Doppler  broadening.  The  resolution  in  this  experiment  is  limited 
by  Doppler  broadening 


Precision  of  linecenter  for  B-X  transition 

doppler  (^BX'°)  Rule  of  thumb  is  that  line  position  can  be  determined  to 

20  within  1/20  the  width  of  the  line  [28:1] 

8a  =  0.004*cm'  * 


To  detect  a  line  position  shift,  the  separation  Aqo, ,  between  two  separately  measured  peaks  must 
be  greater  than  the  uncertainty  of  the  line  position. 

Aa  tot  -=50 

In  wavenumbers 

Aa  tot  “  0.004‘cm"* 
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tot 


Av  =  1.269*10' 


,8  1 


In  frequency 


The  minimum  Doppler  shift  that  can  be  detected  is  then  given  by: 


^BOXO 


Shift  =  1.657*10 


This  is  of  comparable  magnitude  to  the  minimum  Doppler  shift  values  measured  in  this  experiment 
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\  13.  ABSTRACT  (Maximum  200  words) 

! 

I  Recently,  Hong  and  Miller'  used  Fourier  Transform  Emission  (FTE)  spectroscopy  to  determine  the  N2^  mobility 
in  helium  by  relating  the  observed  Doppler  sliift  of  the  N2^  (B-X)  transition  to  the  measured  electric  field  in  the 
discharge  as  determined  b>'  probes.  In  order  to  establish  the  technique  in-house,  this  effort  duplicated  that 
experiment  and  investigated  the  sensitivity  and  utility  of  the  method  Doppler  shift  measurements  in  this 
experiment  disagreed  significantly  with  those  reported  Hong  and  Miller.  Attempts  to  reduce  the  data  (and 
thus  determine  the  mobility)  using  a  model  developed  Hong  and  Miller  were  unsuccessful.  The  model  did  not 
qualitatively  or  quantitatively  describe  the  Doppler  shifts  observed  in  this  work.  At  this  time,  an  explanation  of 
tlie  differences  between  the  two  experiments  is  not  known.  Possibilities  include  effects  of  striations  on  the 
j  Doppler  shift  measurements,  competition  between  one-  and  two-step  excitation  of  the  N2^  (B-X)  state,  and  failure 
I  to  interpret  the  Hong  and  Miller  model  correctly.  Subsequent  attemfits  to  model  Hong  and  Miller’s  own  data 
I  (taken  from  their  figures)  yielded  mobility  parameters  which  disagreed  with  their  published  results. 
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